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ABSTRACT 

Stellar population synthesis (SPS) provides the link between the stellar and dust content of galaxies and their 
observed spectral energy distributions. In the present work we perform a comprehensive calibration of our 
own flexible SPS (FSPS) model against a suite of data. These data include ultraviolet, optical, and near-IR 
photometry, surface brightness fluctuations, and integrated spectra of star clusters in the Magellanic Clouds 
(MCs), M87, M31, and the Milky Way (MW), and photometry and spectral indices of both quiescent and post- 
starburst galaxies at z ~ 0. Several public SPS models are intercompared, including the models of Bruzual 
& Chariot (BC03), Maraston (M05) and FSPS. The relative strengths and weaknesses of these models are 
evaluated, with the following conclusions: 1) The FSPS and BC03 models compare favorably with MC data 
at all ages, whereas M05 colors are too red and the age-dependence is incorrect; 2) All models yield similar 
optical and near-IR colors for old metal-poor systems, and yet they all provide poor fits to the integrated J—K 
and V— K colors of both MW and M31 star clusters; 3) FSPS is able to fit all of the ultraviolet data because 
both the post-AGB and horizontal branch evolutionary phases are handled flexibly, while the BC03 and M05 
models fail in the far-UV, and both far and near-UV, respectively; 4) All models predict ugr colors too red, 
D„4000 strengths too strong and HSa strengths too weak compared to massive red sequence galaxies, under 
the assumption that such galaxies are composed solely of old metal-rich stars; 5) FSPS and, to a lesser extent, 
BC03 can reproduce the optical and near-IR colors of post-starburst galaxies, while M05 cannot. Reasons for 
these discrepancies are explored. The failure at predicting the ugr colors, D„4000, and H5a strengths can be 
explained by some combination of a minority population of metal-poor stars, young stars, blue straggler and/or 
blue horizontal branch stars, but not by appealing to inadequacies in either theoretical stellar atmospheres or 
canonical evolutionary phases (e.g., the main sequence turn-off). The different model predictions in the near- 
IR for intermediate age systems are due to different treatments of the TP-AGB stellar evolutionary phase. We 
emphasize that due to a lack of calibrating star cluster data in regions of the metallicity-age plane relevant for 
galaxies, all of these models continue to suffer from serious uncertainties that are difficult to quantify. 
Subject headings: stars: evolution — galaxies: stellar content — galaxies: evolution 



1. INTRODUCTION 

The spectral energy distribution (SED) of a galaxy contains 
a wealth of information regarding its star formation history, 
dust content, and chemical abundance pattern. These proper- 
ties provide essential clues to the physical processes govern- 
ing the formation and evolution of galaxies from high redshift 
to the present. It is therefore highly desirable to have a robust 
method for extracting the physical properties of galaxies from 
their SEDs. 

The process of translating observed SEDs into physical 
properties is, unfortunately, very challenging because it re- 
quires 1) an accurate understanding of all phases of stellar 
evolution, 2) well-calibrated stellar spectral libraries for con- 
verting stellar evolution calculations into measurable fluxes, 
3) an initial mass function (IMF), specifying the weight given 
to each stellar mass, 4) detailed knowledge of the star-dust 
geometry in conjunction with an appropriate extinction curve; 
i.e., knowledge of the physical conditions of the interstellar 
medium (ISM). Each of these requirements depend on chem- 
ical composition, further compounding the problem. Com- 
bining these ingredients in order to predict the spectrum of 
a galaxy is known as stellar p opulation synthesis (SPS), and 
has a n extensive histor y (e.g., [T inslev & Gunnlll976t IT inslev 
1980; iBruzuall [l98l iRenzini & Buzzonil 119861: iBuzzonl 
1989 f iBruzual & Chariot! 11993b IWorthevI Il994t IMarastonl 
19981 iLeifherer et al.lll999HFioc & Rocca-Vo lmerange ll997r 



I VazdekisI T l999: Yi etai] l200l IBruzual & Charlotl [2003 : 
Jimenez et all 120041: iLe Borgne et alJ l2004t IMarastonl 12005 : 
Schiavonl l2007HCoelho et all 120071: IConroy et alJ 12009a: 
Molla et al l20MlKotulla et al.ll2009h . 

In the past two decades enormous progress has been made 
on each of the requirements mentioned above. Yet, substantial 
uncertainties remain. A non-exhaustive list includes 1) the 
treatment of the core convective boundary in main sequence 
stars (i.e., convective core overshooting); 2) metallicity- 
dependent mass-loss along the red giant branch (RGB) and, 
relatedly, the morphology of the horizontal branch (HB); 
3) the treatment of the thermally-pulsating asymptotic gi- 
ant branch (TP-AGB) phase; 4) the spectral libraries, espe- 
cially for M giant s, TP-AGB stars, and st ars at non-solar 
metallicities (e.g., iMartins & Coelhol |2007[); 4) blue strag- 
gler (BS) stars and th eir ubiquity (e.g., iPreston & S neden 
I2000t iLi &Hanl 120081) : 5) the effects of binarism, which 
might be especially r elevant for massive star evolution (e.g., 
Eldrid ge et al.l 120081) and BS stars; 6) the importance of 
rotati on on massive star evolution (e.g., iMevnet & M aeder 
2000); 7) non-solar abundance ratios, which effect not only 
the stella r evolution calcula tions but also the stellar spec- 
tra (e.g., ICoelho et al.l 120071) : 8) the unknown dependence 
of the IMF on ISM properties such as metallicity and pres- 
sure. These uncertainties can in many cases dramatically im- 
pact the ability to conver t observables into physical proper- 
ties and vice-versa (e.g.. lCharlot"eTaIi ri996: Ch ariol[l996l 
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Yll200l iGallart et al.ll200lt iMaraston et al J 120061 iLee et al. 
20071: IConrov et al.ll2009at iMuzzin et al.ll2009b IConrov et al. 
2009M - 

Thankfully, there exists a wide array of data capable of con- 
straining SPS models. By far the most common type of ob- 
ject used for comparison is the star cluster 1 . The approxi- 
mately uniform age and metallicity of the stars within star 
clusters and the lack of internal reddening (except perhaps 
for very young clusters) affords direct comparison between 
them and the most basic ingredients in SPS — the stellar 
evolution calculations and stellar spectral libraries. Com- 
parisons to star clusters are either made in CMD space or 
in an integrated sense (i.e., the light from all the cluster 
stars are added together). The former technique provides 
a sensitive probe of the main-sequence including the turn- 
off point, sub-giant and red giant branches (e.g., IWortheyl 
[19941: IBruzual & Charlojl2003MAn et al.ll2009l while the lat- 
ter technique provides a more robust measure of the rarer 
brightest stars that dominate the integrated light from the clus- 
ter. It is the latter method of comparison that is more rele- 
vant if the goal of SPS modeling is to understand the inte- 
grated light from galaxies. For this reason, the integrated col- 
ors, surface brightness fluctuations (SBFs), and spectra of star 
clusters have been used extensively to compare and calibrate 
SPS models (e.g.. IBruzual & Charlol 120031: IGonzalez et al I 
2004 |Marastori2q05t ICohen et alj|2007t ICordier et al l | 2007[ 
Pessev et al.ri2008t iMarigo et all 12001 iKoleva et alj l2008t 
Lee et al.ll2009bl 

There are two significant sources of concern when attempt- 
ing to calibrate SPS models with star cluster data. The first 
concern is that the brightest stars, which dominate the in- 
tegrated light, are rare, and thus stochastic effects must be 
carefully modeled. Recently, this issue has been addressed 
observationally by stacking clusters in bins of age in order 
to synthesize 'sup erclusters' that are relatively unaffected by 
stochastic effects dPessev et al.ll2006l IGonzalez et al.l l2004). 

The second concern is less tractable, even in principle, and 
arises from the fact that the primary sources of star clusters for 
which ages and metallicities can be reliably estimated are the 
Milky Way (MW), M31, and the Large and Small Magellanic 
Clouds (LMC and SMC, respectively). This fact imposes se- 
vere restrictions on the region of the age-metallicity param- 
eter space that can be constrained with star cluster data. For 
example, one of the most important regions of this parameter 
space for studying galaxies — old and metal rich — contains 
very few star clusters (except in the bulge of the MW, although 
high extinction diminishes the utility of star clusters there). In 
Local Group star clusters, a bin in age contains clusters of a 
relatively narrow range in metallicity owing to the simple fact 
that the metallicity of galaxies tends to increase with age. 

Because of these issues, entire galaxies are often also used 
to assess the accuracy of SPS models (e.g., IWorthevlll994t 
Bruzual & CharlotH2003T: IMaraston et alJl2006t lEminian et al.l 
2008). Galaxies are not subject to the two concerns mentioned 
above, but using them to constrain SPS models can be very 
difficult because 1) galaxies contain stars of a range of ages 
and metallicities; and 2) starlight in galaxies is attenuated by 
interstellar dust. Nevertheless, if subsamples of galaxies are 
carefully constructed with known and regular properties, then 
galaxies can be used to assess the reliability of SPS models. 

The present work seeks to provide a comprehensive com- 
parison between SPS models and observed star clusters and 

1 Herein both open and globular clusters will be referred to as star clusters. 



galaxies. Such a comparison is timely because of an abun- 
dance of new, high-quality data, including near-IR photome- 
try and SBFs of star clusters in the LMC and SMC, UV pho- 
tometry of clusters in M31 and M87, CMDs and integrated 
spectra of clusters in the MW, and optical and near-IR pho- 
tometry and spectral indices for large samples of low-redshift 
galaxies. We will c onsider not only our own SPS model 
(IConrov et al. 2009a) but also the com monly used mo dels of 
IBruzual & Chariot! (120031 BC03) and lMarastonl(l200l M05). 
Such a detailed comparison between these popular models has 
not been undertaken until now, and will thus provide a much- 
needed evaluation of the relative strengths and weaknesses 
of these models. In addition, the flexible nature of our own 
model will be exploited to explore the extent to which vari- 
ous uncertain phases of stellar evolution, including thermally- 
pulsating asymptotic giant branch stars (TP-AGB), blue HB 
stars, and post-AGB stars, can be constrained by the data. 
This task is critical if we are to have confidence in the derived 
physical properties of galaxies, such as stellar masses and star 
formation rates. 

We proceed as follows. £|2]contains a detailed description of 
the star cluster and galaxy data used to constrain the models, 
which are themselves described in !j3] In fj4]we present an 
extensive comparison between models and data. A discussion 
and summary are provided in |J5] and Sj6] respectively. The 
zero points of the star cluster UBVRIJHK magnitudes are in 
the Veg a system; all other m agnitudes are quoted in the AB 
system (lOke & Gunnlll983l) . Where necessary, we adopt a 
flat ACDM cosmology with (fl m ,Cl A ,h) = (0.26, 0.74,0.72). 

2. DATA COMPILATION 

This section describes in deta il the sources a nd tre atment of 
data both for star clusters ($2J} and galaxies (i )2.2l i. 

2.1. Star cluster data 

The calibration of SPS models against data is essential be- 
cause of the many uncertain model ingredients, including late 
stages of stellar evolution and stellar spectral libraries. As 
mentioned in the introduction, the star cluster is the ideal ob- 
servational datum for SPS calibration because it is, at least 
approximately, a coeval set of stars at a single metallicity. In 
order to provide comprehensive constraints on SPS models, 
one would like to have data on star clusters spanning a wide 
range in age, metallicity, and wavelength. Unfortunately, ow- 
ing to the details of the chemical history of the MW and its 
satellites, the assembly of such a dataset from nearby galax- 
ies, including our own, is challenging. 

In the present work attention is focused on data from the 
Milky Way (MW), M31, and the Magellanic Clouds (MCs). 
The first two galaxies are sources primarily of old, metal-poor 
clusters, while the MCs provide clusters spanning a wider 
range of ages and metallicities. 

Before describing the datasets in detail, it is worth em- 
phasizing those areas of parameter space that are not well- 
sampled and yet are important for interpreting the integrated 
light from galaxies. Massive star clusters (M > 10 4 M Q ) at 
solar or super-solar metallicities are not common in the Lo- 
cal Group. Less massive star clusters of such metallicities are 
more common, but stochastic effects complicates model com- 
parison. The star clusters NGC 6791, NGC 188, and the bulge 
star clusters NGC 6553 and NGC 6528 are notable exceptions 
in that they are relatively old and metal-rich. These clusters 
will be used herein for calibration. 



SPS Model Calibration, Comparison, and Evaluation 



3 



2.1.1. The Milky Way 

Data on MW star clu sters ar e taken primarily from the glob- 
ular cluster catalog of iHarrisI (119961) . This catalog provides 
UBVRI photometry, [Fe/H] metallicities, and E(B-V) red- 
dening values for 150 clusters. The photometry is corrected 
for Galactic extinction using the re ddening values prov ided 
in the catalog in conj unction with a ICardelli et al.l dl989l) ex- 
tinction curve. The Harris (1996) catalog is complemented 
with JHK photometry from th e Two Micron All Sky Survey 
(2 MASS:lSkrutskie et"aT1l2006l) for 106 clusters as compiled 
bv lCohen et al.l (120071) . For most of the comparisons between 
models and data the MW sample is restricted to clusters with 
E(B-V) < 0.2 in order to reduce the uncertainties associated 
with de-reddening. When presenting V — K and J—K colors, 
all clusters with E(B-V) < 0.4 are included since reddening 
c orrections are smaller in the near-IR. 

Schi avon et al.l (120051) have measured integrated, flux- 
calibrated, high signal-to-noise spectra of 40 MW star clus- 
ters spanning a range of metallicities and HB morphologies. 
The spectra span the wavelength range 3360 < A < 6430A 
at a FWHM resolution of 3. 1 A and samplin g of 1 .OA . These 
spectra have been de-reddened using the lCardeUietal] fl989) 
extinction curve with E(B-V) values tabulated by Schiavon 
et al. Metallicities are also provided by Schiavon et al. and 
will be used herein, except for 19 clusters that have updated 
[Fe/H] measurements from lCarretta et all (I2009I) . 

2.1.2. M31 

Photometry of star clusters in M3 1 is available in th e 
Revised Bologna Catalog (RBC v3.5; iGalleti et all 120041) . 
Available photometry includes UBVRI from a variety of 
sources, JHK from 2MASS, and new ultraviolet observations 
from the Galaxy E volution Explore r {GALEX; iMartin et al.1 
120051) presented in iRey et all (120071) . Spectroscopic metal- 
licities of a subset of M31 star clusters are taken from 
iBarmbv et all d2000l) and lPerrett etaTI (120021). Galactic E(B- 
V) reddening estimates are derived from lSch legel etaT1 cl998) 
via t he utilities in the kcorr ect v4 . 1 . 4 software pack- 
age dBlanton & Rowe is 2007). Photometry is corrected for 
Galactic extinction v ia E(B - V) and the extinction curve of 
ICardelli et al.l (119891) . In addition to Galactic reddening, the 
colors of M31 star clusters will also be reddened by dust 
within M3 1 . Unfortunately it is not straightforward to make 
such a correction, and therefore no additional correction will 
be applied. In later sections attention will be focused on the 
near-IR colors of M31 clusters, where dust corrections should 
be small. 

The RBC contains a number of quality flags. Clusters are 
selected to be definite globular clusters (classification flag 
=1), with a definite confirmation (confirmati on flag =1), and 
to not be in the young cluster catalog of Fus i Pecci et al.l 
(2005). This sample is further restricted to clusters with 
E(B-V) < 0.2 to reduce the uncertainties associated with 
Galactic reddening. 

The final cluster catalog, which includes the aforemen- 
tioned quality cuts and metallicity measurements, contains 22 
clusters, 20 of which include GALEX NUV photometry, and 
12 of which include GALEX FUV photometry. 

2.1.3. The Magellanic Clouds 

Data on star clusters in th e MCs come principally from two 
sources. iPessev et al.l d2006l) presented a carefully constructed 
catalog of 2MASS JHK photometry for a sample of 75 star 



clusters. In IPessev et al.l (120081) this near-IR photometry is 
combined with BV photometry from the literature. These au- 
thors construct a 'test sample' of star clusters where high qual- 
ity data are available to estimate main sequence turn-off ages 
from CMDs and metallicities from either the CMDs or indi- 
vidual stellar spectroscopy. These high-quality data are then 
binned by age to produce average colors and metallicities as a 
function of cluster age. 

In the youngest age bin a le ss direct method is used for 
determining cluster ages in the IPessev et al.l (l2008|) sa mple. 
The 'S-parameter', introduced by Elso n & Falll (119851) . pro- 
vides a simple empirical relation between cluster age and a 
combination of integrated U-B and B - V color s. High qual- 
ity age estimates fro m HST observations from iKerber et al.l 
(2007) were used by IPessev et al.l (120081) to update the S- 
parameter-age calibration in order to provide age estimates 
for the youngest clusters in their sample. The Pessev et al. cal- 
ibration is substantially diffe rent from the commonly used re- 
lation of Girardi et al] d 1995b . It is important to keep in mind 
that since ages are based on CMD fitting, they will depend 
on the stellar evolution calculations used to convert main se- 
quence turn-off points into ages. In particular, the treatment 
of convective overshooting has a strong impact on the age- 
main sequence turn-off point relation for young and interme- 
diate ages (t < 2 x 10 9 yrs). This sensitivity arises because 
overshooting increases the amount of fuel available for core 
hydrogen burning. We will return to these issues when these 
data are compared to models. 

The secon d significant source of data on MC star clusters 
comes from Gonzal ez et al.l ([2004). These authors analyze 
2MASS data to produce not only integrated near-IR colors 
but also surface brightness fluctuations in the 2MASS JHK 
filters. Clusters a re stacked accordin g to the classification 
system devised by ISearle et"afl dl980l SWB), which closely 
tracks the S-parameter described above. The re lation betwee n 
SWB type, age and metallicity is adopted from ICohenl (1 1 9 8 2 ) . 
We have independently verified that these ages agree with the 
ages inferred from the S-parameter calibration discussed in 
IPessev et al.l (|2008[). This sample has been augmented with 
V -I colors (Gonzalez-Lopez lira et al.ll2.Q05l), where / m agni- 
tudes are from the DENIS s urvey ( Epchtein et al] 1 19971) and 
V magnitudes are taken from lvan de n Bergh ( 1981). We point 
out that the Pessev et al. and Gonzalez et al. samples are 
not disjoint. Comparison of these two catalogs will provide a 
valuable cross-check on the reliability of the data. 

The novel feature of these two datasets is that the clus- 
ter data have been stacked in bins of age in order to reduce 
stochastic effects. A number of authors have discussed the 
importance of stocha stic effects arising from small numbers 
of ve r y bright stars (lLancon & Mouhcine 2000; IBruzual A.I 
I2002t ICervino & Luridianal 12004, iBruzual & Chariot! 120031). 
Both Lancon & Mouhcine (2000) and Cervino & Luridiana 
(2004) provide limits on the minimum cluster mass that 
should be satisfied so t hat stochastic fluctuations in integrated 
colors are minimized. IPessev et al.l d2008l) demonstrate that 
the masses of their stacked clusters easily satisfy both con- 
straints, ensuring that their stacked results can be robustly 
compared to SPS models. The m asses of the stacked clus- 
ters in t he iGonzalez et al] d2004l) sample are comparable to 
those in IPessev et al.l (120081) . and so their results should also 
be relatively unaffected by stochastic effects. 

Finally, we include for c ompleteness the V- K and J-K 
colors of LMC clusters from lPersson et al.l (fl983). which have 
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been used in many comparisons between observations and 
SPS models. These J and K magnitudes are converted from 
the California Institute of Technology (CIT) photometric sys- 
te m to the 2MASS sy stem using the transformations provided 
by Carpenter (2001). Ages for these clusters were assign ed 
b ased on the S-param eter calibration of Pes sev et al.l ((2008). 

Pess ev et"ail (|2006) presented a detai led compari s on be - 
tween their near-IR colors and those of Pers son et al.l dl983l) . 
in some cases noting substantial differences in the sense that 
colors from Persson et al. were redder than the colors from 
Pessev et al. They attributed this difference to the fact that 
Persson et al. centered their photometry on the brightest part 
of the cluster in order to maximize the measured flux, while 
Pessev et al. center their apertures using a Gaussian smooth- 
ing technique. In some cases the Persson et al. centroids are 
significantly offset from the true dynamical center of the clus- 
ter because of off-center IR-bright sources within the cluster. 
For this reason, in addition to the fact that the data are stacked 
to reduce stochastic effects, the results of Pessev et al. should 
be preferred to the older results of Persson et al. We will re- 
turn to this issue in £j4] 

2.2. Galaxy data 
2.2 A. The utility of galaxies 

For the reasons mentioned in the previous section, star clus- 
ters are the ideal observational data to calibrate SPS mod- 
els. There are however at least two reasons why it is useful 
to also consider SPS constraints from galaxies: 1) galaxies 
tend to be of higher metallicity than old star clusters, and so 
SPS model predictions for old ~ Z Q populations can be as- 
sessed with galaxies composed of old stars; 2) galaxies con- 
tain many more stars than star clusters, and so stochastic ef- 
fects are never a source of concern when comparing models 
to galaxy data. 

These two benefits of using galaxies to calibrate SPS mod- 
els are exploited herein by considering two classes of galax- 
ies. The first class consists of bright red sequence, 'quies- 
cent' galaxies, which contain old metal-rich stars, and will 
therefore be useful for calibrating old metal-rich SPS model 
predictions. The second class are called 'K+A', or 'post- 
starburst' galaxies, as their spectra are dominated by a mix- 
ture of A and K t ype stars and show no signs of ongoin g 
star formation (e.g jDressler & Gunnll98llGoto et al.l2003l) . 
Such galaxies are presumed to have had their star formation 
abruptly truncated. The presence of A stars in their spectra 
indicates that the gala xy contains a significan t population of 
< 1 Gyr old stars (e.g. jLeonard i & Rosa ll996l) . These galax- 
ies will provide valuable constraints on the importance of TP- 
AGB stars, as th e prominence of this stellar phase pe aks for 
ages w 1 -2 Gyr (lMarastonl2005l:IConroy et al.ll2009al) . Since 
stars spend a small amount of time in this phase, observational 
constraints from star clusters can be difficult to interpret be- 
cause of stochastic variations. As mentioned above, a stack- 
ing technique is employed when comparing to star clusters so 
that stochastic effects are minimized. Nonetheless, it will be 
useful to derive independent constraints from galaxies where 
stochastic effects are known to be negligible. 

As mentioned in the introduction, the use of galaxies to as- 
sess the accuracy of SPS models is complicated by at least 
two effects: the attenuation of starlight by dust and the fact 
that galaxies contain stars of a range of ages and metallicities. 
These complications are addressed where relevant in 2] 



2.2.2. Data sources 

The source of all galaxy positions, optical photo metry, and 
redsh ifts is the Sloan Digital Sky Survey (SPSS: lYork et all 
2000), as made available throug h the hybrid NYU Value 
Added Galaxy Catalog 2 (VAGC: iBlaiiton et alJl200l . The 
SDSS photometry is k-corrected and corrected for Galac- 
tic extinction via the kco rrect v4 . 1 . 4 software package 
dBlanton & Row eis 2007). The photometry is k-corrected to 
z = 0.1 - the median redshift of the sample - in order to min- 
imize the magnitude of the k-correction. Spectroscopic in- 
dices are provided by the MPA/JHU database 3 . 

Near-IR photometry in YJHK filters is provided by the 
UKIRT Infrared Dee p Sky Survey 4 (UKIDS S). The UKIDSS 
project is defined in Lawrence et a l. (2007). UKIDS S use s 
the UKIRT Wide Field Camera rWFCA M-lcksali et al.l2 007). 
The photometric system is describ ed in iHewett et al.l (1 2006). 
and the calibration is described in Hod gkin et al.1 (120091) . The 
pipeline processing and scien ce archive are desc ribed in Ir- 
win et al (2009, in prep) and Ham blv et al.l (|2008). We make 
use of Data Release 3, which is publically available. UKIDSS 
photometry is matched to the SDSS with a search radius of 
5". Petrosian magnitudes are used for the UKIDSS photome- 
try. We do not consider colors that are combinations of SDSS 
and UKIDSS photometry, so it is not essential that the photo- 
metric apertures be precisely matched. 

The UKIDSS survey goes much deeper than 2MASS, and is 
therefore more desirable. For example, typical 2MASS pho- 
tometric errors for galaxies at z ~ 0.1 are 0.2 mags, compared 
to ~ 0.01 for UKIDSS photometry. Galacti c extinction is 
corre cted using E(B-V) red dening estimates (ISchlegel et al.l 
1998) in conjunction with a Card elli et al.l (|1989) extinction 
curve. K-corrections for the UKIDSS photometry are esti- 
mated from the kcorrect v4 . 1 . 4 software package us- 
ing the 2MASS JHK filters since those, and not UKIDSS fil- 
ters, are available in the code. The filter transmission curves 
for 2MASS and UKIDSS are very similar and so this approx- 
imation is justified. 

Three galaxy samples are considered . The first are K+A 
galaxies as identified by Oui ntero et al.l ( 120041) . This sample 
contains 1194 galaxies. Of these, only 186 have UKIDSS 
photometry. When considering near-IR colors of the K+A 
sample we use only these 186 galaxies, while when SDSS 
colors are considered we use the full sample. 

The second sample consists of bright red sequence ('qui- 
escent') galaxies. This sample is defined by Mj? 1 < -20.5, 
where the superscript indicates that the magnitude has been 
k-corrected to z = 0.1, and 

(g-r) QA > -0.03 (M°/ l + 23.0) + 0.94, (1) 

which divides galaxies by the valley visible in the color- 
magnitude diagram and effectively isolates red sequence 
galaxies. Applying these cuts and requiring that galaxies have 
both SDSS and UKIDSS photometry results in a sample of 
5340 galaxies at 0.08 < z < 0. 12. 

The third sample is derived from the spectroscopic indices 
catalog. Galaxies are selected with 0.05 < z < 0.1 and with 
a median spectroscopic signal-to-noise of S/N > 25. From 
this catalog two classes of galaxies are defined: those with 
velocity dispersions measured from the stellar continuum of 

2 http : / /sdss . physics . nyu . edu/vagc/ 

3 http : / /www . mpa-garching . mpg . de/SDSS/ 

4 http : / /www .ukidss.org/ 
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100 < a < 200 kms" 1 and those with a > 200 kms" 1 . These 
samples contain 13214 and 6774 galaxies, respectively. 

3. SPS MODEL CONSTRUCTION 

3.1. Flexible SPS 
3.1.1. Theoretical isochrones 

SPS mod el con struction closely follows that of 
IConrov et alj (l2009al) . to which the reader is referred 
for details. In brief, this SPS code (referred to below as 
FSPS, for 'flexible' SPS) combines stellar evolution cal- 
culations with stellar spectral libraries to produce simple 
stellar populations (SSPs), which describe the evolution 
in time of the spectrum of a mono-metallic, coeval set of 
stars. The initial stellar mass function (IMF) speci fies the 
weigh ts given to stars of various masses. The IMF of Kroupa 
(2001) is adopted throughout. FSPS also contains a number 
of phenomenological dust attenuation prescriptions and 
produces spectra for composite stellar populations. Our SPS 
code is open-source and publicly available 5 . 

The novel feature of this SPS model is that the essential 
aspects of SPS are integrated in a flexible way so that the user 
may choose, for example, their preferred set of isochrones and 
stellar spectral libraries. In addition, the user may alter the 
weights given to and physical parameters of various advanced 
stages of stellar evolution (see below), and may specify the 
IMF of their choosing. The details of these various model 
choices are now described. 

We make use of bot h the Padova 6 dGirardi et all l2000t 
IMarigo & Girardl 120071: iMarigo et al.1 120081) and BaSTI 7 
dPietrinferni et al. , 2004; Cordier et al. 2007) stellar evolution 
calculations, both in the versions available as of November 
2009. These two groups follow stellar evolution from the 
main sequence through the thermally-pulsating asymptotic 
giant branch (TP-AGB) phase, although the treatment of the 
TP-AGB phase is very different between the two groups. 

The Padova calculations exist for metallicities in the range 

^10.2 



< t < 10 1U Z yrs, and 
BaSTI outputs cover 



2 x 10" 4 < Z < 0.030, for ages 10 6 ' 6 
for initial masses 0.15 < M < 67 M Q . 
a somewhat narrow range in ages, 10 75 < t < 10 10 2 yrs, and 
initial masses 0.5 < M < 10.0M Q but include higher metal- 
licities: 3 x 10" 4 < Z < 0.040. The BaSTI database pro- 
vides models for both solar-scaled and a-enhanced chemical 
compositions, for two assumptions regarding convective core 
overshooting, and for two stellar mass-loss rates. We use the 
solar-scaled, no overshooting, and r\ = 0.4 mass-loss rate pa- 
rameter model. For both lower and higher masses the Padova 
calculations are stitched onto the BaSTI grid. The Padova 
and BaSTI models are compared in Figure Q] In this figure 
isochrones are compared at several ages and metallicities, and 
in addition the location of the TP-AGB onset is shown. 

In order to follow stellar evolution from the end of the TP- 
AGB phase through the post-AGB phase, the above evolu- 
tionary tra cks have been supplemente d with the post-AGB 
models of Vassiliadis & Wood (1994). The contribution of 
post-AGB stars to the integrated SED is highly uncertain for 
a number of reasons, including a lack of calibrating obser- 
vations and the unknown duration of heavy enshroudment as 
post-AGB stars leave their M giant phase. Recent observa- 
tions of M32 have found that post-AGB stars are much less 

5 www .astro. princeton. edu/^cconroy /SPS/ 
6 http :// stev . oapd . inaf . it /cgi-bin/cmd 
7 http : / /www . oa-teramo . inaf . it /BASTI 
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FIG. 1 . — Theoretical HR diagram comparing the BaSTI and Padova evolu- 
tionary calculations at several ages and metallicities. Isochrones are plotted 
from the main sequence through the TP-AGB phase. The onset of the TP- 
AGB phase is marked by solid symbols along each isochrone. The 10 10 yr 
solar metallicity BaSTI isochrone and the 10 8 yr Padova isochrones do not 
contain TP-AGB stars. The difference between the models at young ages 
is due to the different treatments of convection in the two models (i.e., the 
handling of convective core overshooting). 



frequent than standard models predict dBrown et al.l 120081) . 
Nonetheless, the inclusion of post-AGB tracks may be im- 
portant for th e interpretation of the ultr aviolet spectra of old 
systems (e.g..Greg gio & Renzinil fl990). 

The FSPS code contains a number of parameters that con- 
trol the weights given to various uncertain phases of stellar 
evolution. For example, the bolometric luminosity and ef- 
fective temperature of TP-AGB stars can be arbitrarily mod- 
ified. These modifications are applied as overall shifts in 
log(Lboi) and log(r e ff) for the entire TP-AGB phase, and are 
pa rameterized by the v ariables A^, and At (see discussion 
in IConrov et al.l l2~009ah . In addition, an extended blue HB 
can be included 8 , as can a population of BS stars. A blue 
HB population is parameterized as the fraction, /bhb of HB 
stars that are in a blue/extended component. Such a flexi- 
ble treatment is motivated by the great variety of observed 
HB morphologies in MW clusters, even at a given metallicity 
dPiotto et al. 2002). These stars are uniformly populated from 
the red clump to log(T e ff/K) = 4.2. BS stars are uniformly 
populated from one half magnitude above the main sequence 
turn-off to two and a half magnitudes br ighter than the turn - 
off, along the zero-age main sequence dXin & De ng 2005). 
The weight given to BS stars is specified as the number of BS 
stars per unit HB star, Sbs ■ The weight given to the post-AGB 
phase can also be arbitrarily modified. The effect of these un- 
certain aspects of stellar evolution will be investigated in later 
sections. 

3.1.2. Spectral libraries 

Stellar spectral libraries are required to convert the out- 
puts of stellar evolution calculations — surface gravities, g, 
effective temperatures, T e ff, and metallicities, Z, — into ob- 
servable spectral energy distributions. Two spectral libraries 
are considered herein. The first is the empirically-ca li brated 
theoretical BaSeL3 .1 library dLejeune et alJ 1 1 9971 11998b 
iWestera et al.l 120 02). supplemented with empirical, average 

8 When we consider blue/extended HB populations we do not add new 
stars but rather move stars from the red clump HB to the blue component. 
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TP-AG B spectra from the library of La ncon & Mouhcinel 
(2002). The BaSeL library was constructed from the model 
atmosp heres of Kurucz (1 9 95, private communication to R. 
Buser). iFluks et all (11994 iBessell et alJ ([1989, 199lJ), and 
lAllard & Hauschildd7l995h . In the TP-AGB library, spec- 
tra are binned both by near-IR color and by whether the star 
is carbon-rich or oxygen-rich. Stars hotter than 50,000^ are 
assumed to be pure blackbodies. This combined library cov- 
ers the full wavelength range (ultraviolet through infrared) at 
a sampling of 10A in the ultraviolet and 20A in the optical. 

The second spectral library cons idered herein is derived 
from the empirical Miles library (Sanchez-Blazque z et al.l 
2006). This atlas contains spectra of ss 1000 stars covering the 
wavelength range 3600 < A < 7500A at a sampling of 0.9A 
per pixel (FWH M» 2.3A). De rived stellar parameters are pre- 
sented in lCenarro et al.l(| 2007). Integrating empirical libraries 
into SPS models is a non-trivial task because the coverage in 
the stellar parameters g, T e ff, and Z is very inhomogeneous. In 
fact many regions of this parameter space are not covered by 
empirical libraries at all, especially at sub-solar and super- 
solar metallicities and for spectral types earlier than A. An 
additional problem with the Miles library is that the spectra 
have been arbitrarily normalized to unity at 5350A. This is 
particularly problematic because the wavelength range is not 
broad enough to compare the integral of the flux to the bolo- 
metric luminosity for a star with a given g, T e ff, and Z. There- 
fore, an additional library must be used to provide the abso- 
lute flux scale. The BaSeL library is used for this purpose. 
The BaSeL library is also used to provide the relative trends 
of spectra with stellar parameters for regions of parameter 
space where the relative trends cannot be self-consistently 
determined from Miles. The Miles library is binned and, 
when necessary, linearly interpolated onto the Kurucz/BaSeL 
log(g)-log(r e ff) grid. This interpolated library is then in- 
terfaced with FSPS in a manner identical to the BaSeL li- 
brary. Five metallicity bins are defined so that roughly similar 
numbers of stars occupy each bin. The bin boundaries are: 
log(Z/Z ) = [-2.0,-1.0,-0.5,-0.15,0.15,0.75]. The result- 
ing average metallicities within these bins are log(Z/Z ) = 
[-1.4,-0.7,-0.3,0.0,0.2]. 

At metallic ities less than solar , the Miles library stars are 
a-enhanced dMilone et al.ll2009l) . At [Fe/H]< -1 the library 
is on average enhanced in Mg by +0.4 dex. This will be rele- 
vant in our comparison to the spectra of MW star clusters in 
later sections. 

There are several important caveats to these libraries that 
the reader should keep in mind. In the empirical TP-AGB 
library, the metallicities of the stars are unknown. The au- 
thors advocate assuming that they are all solar metallicity 
since a significant fraction are located in the solar neighbor- 
hood (although many are also located in the MCs). In order 
to use the library at non-solar metallicities, the authors sug- 
gest the use of the theor etical metallicity-d ependent color- 
temperature relations in Bess ell et al.l (Q991) in order to con- 
vert the observed I-K colors into temperatures at other metal- 
li cities. 

IWestera et al.l (120021) found that the original theoretical 
atmosphere calculations used in the BaSeL library could 
not simultaneously match observed UBVRIJHKL color- 
temperature relations for individual stars and CMDs of glob- 
ular clusters. The reason for this was thought to be inade- 
quacies in the stellar evolution calculati ons used in the CMP 
comparisons. Because of this problem, Wes tera et al.1 (120021) 



presented two normalizations of the BaSeL library, one which 
matched the color-temperature relations, and the other which 
matched the globular cluster data. These two normalizations 
are referred to as the 'WLBC' and PADOVA2000' libraries, 
since the Padova evolutionary calculations were used in the 
globular cluster comparisons. The 'WLBC' normalization 
will be used herein because this calibration is not tied to any 
particular isochrone set. 

In summary, we have in hand a flexible SPS model from 
which we can choose two sets of isochrones (Padova or 
BaSTI) and two sets of spectral libraries (BaSeL or Miles). 
In addition, there are a number of free parameters character- 
izing various uncertain aspects of advanced stages of stellar 
evolution. These various options will be exploited in later 
sections. 

3.1.3. Modifications to the Padova isochrones 

For the Padova isochrones, the TP-AGB phase has been 
modified in order to generate acceptable agreement with the 
data. Specifically, the TP-AGB phase has been calibr ated 
again st bo th MC star clusters and post-starburst galaxies (i !4.2l 
and i !4.6l respectively). These calibrations result in an ad- 
justment of the TP-AGB parameters Ai an d A j (quantifying 
overall shifts in log(Lboi) and log(r e ff); see S Q.U as follows: 

( -1.0+(f-8.0)/1.5, 8 < r < 9.1 

L ~\ -max{min{0.4,-log(Z/Z Q )},0.2}, r > 9.1, 

(2) 

a f/^'y/^X- 025 8 <r<9.1 m 
A T = < 1 0.0, otherwise (3) 

[ +0.1-min{(f-9.1)/1.5,0.2}, r > 9.1, 

where t is the age of the population measured in log(yrs). We 
emphasize that the values of Ai and A 7- chosen to fit the data 
are not unique, and that, unless explicitly stated in the above 
definitions, they are independent of metallicity. Throughout, 
reference to the 'modified Padova model' refers to these mod- 
ifications. 

3.2. Additional SPS models 

We now describe the salient aspects of the BC03 and M05 
SPS models. 

The BC 03 model util i zes the Padova isochrones circa 
1994 (lAlongiet al.1 Il993t iBressan et al1l!993t iFagotto et al l 
11994 . supplemented with the IVassiliadis & Wood! (Il993h 
model s for TP-AGB stars and the Vassiliadis & Wood! 
dl994 models for post-AGB star s. BC03 use the models 
of iGroeneweeen & de Jong] (1 1993b to differentiate between 
carbon-rich and oxygen-rich TP-AGB stars. There are no 
blue/extended HB stars in the BC03 models. Their models are 
constructed with both the theoretical BaSeL library (described 
in the previous section ) and the empirical STELIB library 
dLe Bo rgne et al. 2003). Spectra for carbon stars are t aken 
from model atmosphere calculations (Hofn eret al.ll2000l) . 

The empirical STELIB library contains optical spectra for 
249 stars at a resolving power of R « 2000 (FWHM- 3A). 
The effective temperatures in the STELIB library are not re- 
liable and so BC03 adopt the color-temperature scale from 
the BaSeL library in order to relate the stellar evolutionary 
outputs to the empirical library. There are very few stars at 
non-solar metallicities and so the reliability of models con- 
structed with thi s library at non-sol ar metallicities is question- 
able. Indeed, iKoleva et al.l d2008) compared several differ- 
ent SPS models and found that the BC03 model constructed 
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FIG. 2. — SSP colors as a function of age, comparing a variety of SPS models to observations. Data are derived from MC star clusters. The model SSPs include 
the predictions from M05 (dot-dashed lines), BC03 (dot-dot-dot-dashed lines) and two models derived from the present work. Using our FSPS code, a model 
was constructed usi ng the BaSTI (dashed lines) and the modified Padova stellar evolution calculations (solid lines). The MC data are from Persson et al] 119831 
P83), Gonzalez-Lopezlira et al. 12005, G05) and Pessev et al. (2008, P08). For the latter two data sets, each data point represents an average over many observed 
star clusters within each age bin. The Persson data is believed to be less reliable for reasons described in the text. The average metallicity of the data decreases 
as the age increases, from log(Z/Z0) = -0.3 at log(f)=7.5 years to log(Z/Z0) = -1.5 at log(f)=10.0 years. The model predictions reflect this change in average 
metallicity with age. 



with STELIB produces unreliable results at non-solar metal- 
licities. 

We make use of the BC03 model that was constructed with 
the STELIB library for the optical spectra. For both ultra- 
violet and near-IR predictions, the STELIB library has been 
extended with the BaSeL library. 



The M05 model is built from the following ingredients. 
Stellar evolution calculations through the main sequence 
turnoff are adopted from Cassisi et al. (JCassisi et alJll997L 
2000). The fuel consumption theorem (Renzini & Buzzoni 
Il986t iMaraston 119981) is then used to incorporate advanced 
stages of stellar evolution through the TP-AGB phase. The 
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FIG. 3. — SSP colors as a function of age, demonstrating the sensitivity to 
the TP-AGB phase. The Padova stellar evolution calculations are used for 
the model predictions. The modified model shown in Figure^(dashed lines) 
is compared to the unmodified Padova calculations (solid lines). Dotted lines 
correspond to a ±0.2 dex variation in metallicity about the modified model. 
The data are from star clusters in the MCs, as in Figure[2] Notice the different 
y— axis ranges in the bottom panels compared to FigureTI] 



evolutionary calculations are converted to observables with 
the theoretical BaSeL spectral library for all stars except 
TP-AGB stars. For thes e stars the empirical library of 
lLancon & Mouhcinel d2002l) is used (this is the same library 
used in our FSPS code). 

Model magnitudes and spectral indices 9 are computed self- 
consistently from the model spectra. Comparisons between 
models in later sections are therefore not plagued by poten- 
tially different filter transmission curves or index definitions. 

4. SPS MODEL CALIBRATION & COMPARISON 
4.1. Overview 

We now turn t o a com parison between observations and 
SPS models. In tj4-.2444-.3l we focus on star clusters, both in 
the MCs, where constraints are strong for intermediate age, 
sub-solar metallicities, and in the MW and M31, where con- 
straints are str ong for o ld metal-poor systems. We then turn 
to galaxies in i j4.4144.6l Two types of galaxies will be consid- 
ered: massive red sequence galaxies and post-starburst galax- 
ies. The former class will provide a qualitative assessment of 

9 Spectroscopic indices are defined by the wavelength intervals specified 
in the MPA/JHU database. The intervals are similar to the classic Lick/IDS 
system ( Worthev et al. 1994), although no attempt has been made to ensure 
that our indices are on the Lick system scale. 



the accuracy the SPS models in the solar to super-solar metal- 
licity, old age regime. The post-starburst galaxies will pro- 
vide an assessment of the models in the approximately solar 
metallicity, intermediate age regime. As discussed in previous 
sections, it is challenging to use galaxies to assess the accu- 
racy of SPS models because galaxies are composed of stars 
with a range of ages and metallicities and dust reddening can 
be importan t. These complications are discussed extensively 
in 34~44l4~6l 

4.2. Constraints from the Magellanic Clouds 
4.2.1. Colors 

Figure [2] shows the colors of MC star clu sters as a func- 
tion of cluster age. As discussed in £12.1.31 individual clus- 
ters are binned in age to produce 'supercl usters' that are rel 
atively immune to stochastic flu ctuations (IPessev et alj|200l 
Gonzalez-Lope zlira et al.ll2005l) . The data from Persson et al 
(1983) are also included, and have n ot been binned in age. 
These data are of inferior quality (see ^2.1. 3t and are only in- 
cluded for comparison bec ause this data is f requently used to 
constrain SPS models (e.g., Maraston 20051). 

A variety of SPS models are compared to the MC star clus- 
ter data in Figure [2] From our FSPS code we have gen- 
erated SSP color evolution for both the Padova and BaSTI 
isochrones. In thi s figur e we make use of the modified Padova 
isochrones (see i j3.1.3t . The models of BC03 and M05 are 
also included. All model predictions reflect the fact that the 
average metallicity of the data varies with age. That is, we 
have interpolated model predictions so that they vary with 
metallicity in as in the data, where the age-metallicity relation 
in the MC data is log(f/yrs)= (7.0,7.7,8.7,9.0,9.5,9.8, 10.1), 
logfZ/Z^) = (-0.28, -0.28,-0.28, -0.28, -038, - 0.68, -1.5) 
(Pess ev et alJl2008tlGonzalez-L6pezlira et alj|2005l) . 

The M05 models are significantly redder than the data at 
intermediate ages , owing to the diffe rent TP-AGB treatment 
in M05 (see also lPessev et al.ll2008l who reach similar con- 
clusions). The M05 colors also have an age-dependence that 
does not agree with the data. The MC cluster ages were es- 
timated via CMD fitt ing using isochrones that include con - 
vective overshooting ( Ke rber et al.ll2007l iPessev et al.l l2008), 
whereas the M05 models do not include overshooting, and 
so one might expect some mismatch in age between the M05 
models and data. However, the differences induced by the 
treatment of overshooting on the MC ages is < 0.2 dex at 
ages 10 8 < t < 10 9 yrs, whereas the discrepancy in the age- 
dependence between M05 and the data is > 0.5 dex in age. We 
therefore regard most of this discrepancy as a shortcoming of 
the M05 model, with only a fraction being due to different 
definitions of cluster age. 

The Persson et al. data are substantially redder than the 
newer datasets from P essev et al. and Go nzalez-Lopezlira 
et al. As discussed in IPessev et alJ (120061) . this is probably 
due to the fact that the Persson photometry was probably cen- 
tered on the brightest star in the cluster, which are often car- 
bon stars, and in any case are very red, and so the photometry 
will be biased red. The difference between Persson's data and 
the 2MASS-based data highlights the care that must be taken 
to construct datasets that can be meaningfully used for SPS 
model calibration. 

The M05 models were calibrated against color-color plots 
of MC star clusters using the Persson et al. data and very wide 
age bins. Stochastic effects will tend to smear star clusters 
along an underling color-color relation because if a cluster is 
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FIG. 4 . — SSP surface brightness fluctuation (SBF) magnitudes and colors as a function of age. Models are compared to data from star clusters in the MCs 
I Gonzalez-Lopezlira et al. 2005, G05; symbols). The models are based on our FSPS code using the Padova and BaSTI isochrones. For the Padova isochrones, 
predictions are included both for the unmodified and modified isochrones. 



scattered redder in one color, it will be redder in another color 
as well. The M05 model calibration is therefore sensitive to 
stochastic effects, and is in addition plagued by the systemat- 
ics in the Persson data. Moreover, the wide age bins used in 
M05 explains why the discrepancies in age seen for the M05 
model in Figure[2]were not seen in their own calibrations. 

The BC03 model and BaSTI isochrones both fare well 
except for H - K and J — K colors at intermediate ages 
(log(f/yrs) « 9.5), for which models are too blue by ks 0.1 
mags compared to the data. This disagreement is probably 
due to the simplistic treatment of TP-AGB stars in these 
models. The BaSTI isochrones also produce redder colors 
at young ages compared to the Padova isochrones, which is 
probably due a lack of convective core overshooting in the 
BaSTI isochrones adopted herein (see Figure[T|i. The data do 
not definitively favor one model over another at log(f/yrs) < 
8.0. 

The disagreement in V - 1 between the models and data at 
intermediate ages is also noteworthy. The models perform 
better at old ages, although it is unclear if this is because of 
the old ages or the lower metallicities of the old star clusters. 
We will return to issues with V—I colors in later sections. 

FSPS model predictions using the Padova isochrones are 
explored further in Figure [3] In this figure we compare the 
unmodified Padova isochrones to the modified model, and we 
also show the effect of a variation in metallicity of ±0.2 dex 
about the modified model. For ages and colors this variation 



with metallicity does not bracket the modified model because 
some metallicity-color relations in the models are not mono- 
tonic and in fact some colors are reddest at approximately 
the average metallicity of the LMC (log(Z/Z Q ) w -0.3). It 
is clear that use of the unmodified Padova models results in 
substantial disagreement with observed star clusters. This dis- 
agreement is driven to a significant extent by the large popula- 
tion of luminous carbon stars predicted in the Padova calcula- 
tions. In our modifications to the Padova models we have sub- 
stantially lowered the bolometric luminosity along the TP- 
AGB phase, which has the effect of down-weighting the im- 
portance of all TP-AGB stars, including carbon stars. 

It is worth discussing in more detail the substantial dis- 
agreement between the unmodified Padova calculations and 
the star cluster data. The Padova group p rovided comparisons 
between their model and observations in Marigo et alj (12008b 
for V — K and J—K colors, finding acceptable agreement. The 
explanation for this difference is twofold. 

First, the data used in their comparison is the old Persson 
data, which suffers from the issues discussed above. Second, 
and perhaps as important, the conversion from L\, Q \ and T e s to 
S EDs was achieved with the theoretical atmosphere models 
of iLoidl et alJ (1200 lh coupled to dust radiative transfer mod- 
els ^from|Groenewegen (2006). Aside from the fact that the 
dust models were never demonstrat ed to fit the V — K colors 
of AGB stars, the ILoidl et al] d2001h carbon star models are 
substantially bluer (by > lmag in V — K) than observed car- 
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FIG. 5. — K— band SBFs as a function of integrated V — I color. Model 
predicti ons (lines) are compared to observations of star clusters in the 
MCs I Gonzalez-Ldpezlira et al. 2005, G05), Elliptical galaxies in the For- 
nax clus ter (Liu et al. 2002, L02), and Ellipticals in a variety of environments 
iJensen et al. 2003, J03). Model predictions are displayed as lines of constant 
age (solid lines; labeled in units of log(f/yrs) in the figure) and lines of con- 
stant metallicity (dashed lines); labeled in units of log(Z/ZQ) in the figure). 
The age-metallicity grids in the top two panels are the same. SSP models 
are constructed with the unmodified Padova isochrones (top panel), modified 
Padova isochrones (middle panel), and the BaSTI isochrones (bottom panel). 



bon stars (R. Gautschy, private communication). We believe 
that our transformation to SEDs is more robust since we use 
observed TP-AGB spectra. Our q ualitative conclusions are 
echoed by Gullieuszi k et al.l (|2008) who compared the popu- 
lation of O-rich and carbon stars in the Leo II dwarf galaxy 
to the Padova models and found that the models predict sub- 
stantially more, and substantially more luminous O-rich stars 
than are seen in that galaxy. Finally, the Padova group uses 
the carbon star luminosity function in the MCs as a constraint 
on their TP-AGB model, but this observable can only provide 
valuable insight if the star formation histories (SFHs) of the 



MCs are known to high precision in the age range relevant for 
carbon star production 10 8,5 - 10 9 15 yrs). 

4.2.2. Surface brightness fluctuations 

Figures |4] and [5] compare near-IR SBFs between several 
FSPS-based models and MC data. In Figure |4] the evolu- 
tion of SBF magnitudes and SBF colors are compared to our 
FSPS model with both the BaSTI and Padova isochrones. We 
show results both for the unmodified and modified Padova 
isochrones. The BaSTI isochrones provide an acceptable fit to 
the /T-band SBFs, but produce SBF colors that are somewhat 
too blue compared to the data. This implies that the AGB and 
TP-AGB temperatures are too hot in the BaSTI calculations. 

The unmodified Padova results show similar discrepancies 
as seen in Figure[3] SBF magnitudes and colors are too bright 
and too red, which can be traced back to the over-abundance 
of carbon stars and the over-lumi nous T P-AGB evolution- 
ary phase. As discussed earlier in £14.2.11 these problems are 
rectified by altering the overall luminosity and temperature 
of stars in the TP-AGB phase. These modifications result in 
SBF magnitudes and colors that are in better accord with the 
data, although /f-band SBFs tend to be somewhat lower than 
the observations at intermediate ages. 

In Figure we compare models and data in /f-band SBFs 
versus integrated V -I colors. In addition to MC star clus- 
ter data we include data from Elliptical galaxies dLiu et al.1 
I2002tl Jensen et al.ll2003l) . It is clear that the modified Padova 
isochrones span the range of data much better than the un- 
modified isochrones. In addition, notice that there is a sub- 
sample of Ellipticals that have SBFs fainter than predicted 
by the models. This has a direct and potentially important 
impact on the inferred mass-to-light ratios of these systems 
and consequently on their inferred stellar masses. SBFs in 
the red/near-IR will correlate inversely with derived stellar 
masses, and therefore if SBFs are overestimated in the mod- 
els (as is the case here), then derived stellar masses will be 
underestimated. 

The models shown in Figure [5] assume solar-scaled 
chemical compositi ons, while many Elliptical galaxies ar e 
q-enhanc ed (e.g.. IWorthev et al.lll99l iThomas et aD 2005). 
However, iLee et al.l d2009al) have shown that the level of 
a-enhancement has almost no effect on the model predic- 
tions in the V -I vs. F160W SBF plane (F160W is similar 
to the //-band), and therefore our conclusions should not be 
effected by these details. 

4.3. Constraints from the MW & M31 

4.3.1. Optical and near-IR colors 

SPS model predictions for the metallicity-dependent opti- 
cal and near-IR colors of old populations are compared to 
observed star clusters in Figure [6] Model predictions are re- 
stricted to ages of 12.5 Gyr, although results are not sensitive 
to this particular age. FSPS models with both the BaSTI and 
modified Padova isochrones are included, as are the predic- 
tions from M05 and BC03. Comparisons to star clusters in 
both the MW and M3 1 are included. For all colors except for 
J-K the MW and M31 data agree (see below), although the 
scatter in M3 1 cluster colors is larger. This is probably due 
to the effects of internal reddening in M3 1 that has not been 
accounted for. 

The models reproduce the trends in the observed colors 
farely well, although there are areas of noteworthy disagree- 
ment. The U — V and B-V colors are somewhat too blue 
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FIG. 6. — SSP colors as a function of metallicity, for old ages (f fs 13 Gyr). Model predictions are compared to observations of globular clusters in the MW and 
M31. The model SSPs include the predictions from M05 (dotted lines), BC03 (dot-dashed lines), and predictions from the present work where either the Padova 
or BaSTI stellar evolution calculations have been assumed. In general the agreement between models and data is acceptable, though the data tend to be redder 
than the models as the metallicity increases. The difference in J— K color between MW and M31 star clusters is striking. 



in the models at log(Z/Z Q ) > -1.0. The FSPS model with 
the BaSTI isochrones perform worst of the models consid- 
ered. We can check whether or not this general disagree- 
ment between models and data in UBV filters is due to the 
adopted spectral library by considering the empirical Miles 
stellar spectral library. Due to the limited wavelength range 
of this library, we can only construct the B-V color. The 



B—V colors from the Miles empirical library are almost indis- 
tinguishable from the colors using the BaSeL library, imply- 
ing that the stellar atmospheres are not the source of the dis- 
crepancy. A similar test was perfor med in BC03 , where they 
compare the BaSeL, STELIB, and lPicklesI (119981) libraries and 
find negligible differences in B-V at Z = Z Q . 
A second possible explanation for the differences between 
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FIG. 7. — CMDs for two metal-rich MW star clusters, NGC 6791 (top rows) and NGC 188 (bottom rows), compared to FSPS model predictions using the 
BaSTI isochrones. Notice the strong BS population above the main sequence turn-off for both NGC 6791 and NGC 188. The adopted cluster parameters are 
shown in the legends. 



models and observations in the U — V and B-V colors is the 
UBV filter transmission curves. The B filter used for these 
data actually consist of two sl ightly different filters (see the 
discussion in iMarastonl 12005). It is difficult to determine 
which filter was used in each observation. The difference be- 
tween these two filters varies from 0.02 to 0.04 mags. The 
U — V color is constructed by adding the U-B and B-V col- 
ors, so this color will also be subject to this source of uncer- 
tainty. The offsets between models and data is of order 0.05 
mags — the same order as the differences induced by the two 
B filters — and so we cannot make stronger comparisons be- 
tween models and data until more reliable optical photometry 
is available for MW and M31 star clusters. 

The morphology of the HB cannot explain the observed dis- 
crepancy because in all SPS models shown in Figure [6] adopt 
a red HB at log(Z/Z ) > -1.0. The HB therefore cannot be 
modified to produce even redder model colors. 

The model V — I colors compare favorably with the data at 
log(Z/Z Q ) > -1.0. At lower metallicities the models diverge 
from the data, with the FSPS model using BaSTI isochrones 
performing worst. 

All models fail at predicting the observed V-K and J-K 
colors, both at low and high metallicities for the latter and 
at higher metallicities for the former. Recall that the near- 
IR data comes fro m a recent analysis of 2M ASS data by 
ICohen et al.l (12007b . Our results echo those of ICohen et alj 
(120071) who compared MW star cluster J-K colors against 
several SPS models. In addition to the M 05 model, they also 
compare to the models of 1Wortheyl(ll994l) and lBuzzonT ( 1989) 



and find similar levels of disagreement with these other mod- 
els; 

ICohen et alj ((2007) speculate that the cause of the discrep- 
ancy could be either 1) inadequate calibration of the stellar li- 
braries or 2) issues with the near-IR filter transmission curves 
(Cohen et al. used published predictions for Johnson JHK 
photometry from these models). Since we compute 2MASS 
magnitudes self-consistently from the model spectra, the sec- 
ond explanation cannot be the source of the discrepancy. The 
problem is not likely to be due to inadequacies in the theoreti- 
cal stellar spectral libraries, since the col ors of giants compar e 
well with observed stars in the near- IR (IWestera et alJ l2002). 
The temperature of the RGB need only be lowered by 200K in 
order to achieve better agreement in V-K (although a lower 
RGB temperature woul d only exacerbate the disagreement in 
J-K). As explored in iH.3.21 one need only change the RGB 
at the brightest ~ 1 mag to in order to produce substantial 
changes in V-K colors. Without definitively identifying the 
source of this problem, we are forced to simply reiterate the 
conclusions of Cohen et al. that near-IR colors of old popu- 
lations cannot be used to constrain the metallicities or ages of 
either GCs or distant galaxies. 

It should also be noted that the near-IR photometry pre- 
sented in Cohen et al., which is based on data from the 
2MASS survey, differs substantially from previous near-IR 
photometry of MW clusters provided by M. Aaronson et al. 
(1977; unpublished). The differences in J-K are 0.1-0.2 
mags in the sense that the latest results are bluer than the older 
data. Since previous models were calibrated against the older 
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FIG. 9. — Spectral indices versus metallicity for old populations. Indices 
are derived from th e integrated spectra of MW star clusters presented in 
ISchiavo n et al. 1 2005, symbols), and are compared to indices of SSPs at t = 10 
Gyr for several SPS models. Here, our FSPS model uses Padova isochrones 
and the empirical Miles stellar library. Additional FSPS predictions that in- 
clude a fraction of blue HB stars (/bhb = 0.5) and BS stars (Sbs - 2) are 
shown. The BaSTI isochrones (not shown), are similar in most respects to 
the Padova isochrones, except for the metal-poor HSa strengths, which are 
larger in the BaSTI models because of a bluer HB. The model predictions of 
BC03 are also included in the figure. 



FIG. 8. — Top panel: Cumulative fraction of V— band light as a function of 
apparent V— band magnitude for a 7 Gyr SSP with [Fe/H]= +0.0. The apparent 
V-band magnitude along the x— axis is scaled by the distance modulus and 
reddening value appropriate for the star cluster NGC 188 shown in FigurefT] 
Red clump HB stars appear at V ~ 12. Bottom panel: Integrated color for 
all stars fainter than V, normalized to the integrated color for all stars. The 
model SSP is the same as in the top panel. Integrated colors in the red and 
near-IR can only be constrained from CMDs if the very brightest and rarest 
stars in the cluster are identified. 



data (e.g., Maraston 2005), it is not surprising that the models 
are all redder than the Cohen et al. data. 

The most dramatic discrepancy in Figure [6] is between the 
J-K colors in MW and M31 star clusters. The M31 clusters 
are approximately 0.2 mags redder in J-K than the MW clus- 
ters. Internal reddening in M31 cannot be the source of this 
discrepancy because reddening effects are quite small in J-K, 
especially in comparison to the colors considered in Figure 
|6l which do not show such an offset. Barmby et al. (2000) 
have also presented JHK photometry of old M31 star clus- 
ters. Their results are consistent with the M31 results from 
the RBC catalog shown in Figure [6j suggesting that the dis- 
crepancy in J- K is not due to data reduction issues. We can- 
not provide a compelling explanation for this discrepancy, but 
note that without an explanation it is impossible to calibrate 
the old metal-poor models to better than ~ 0.2 mags in the 
near-IR. 



4.3.2. Metal-rich CMDs 

Calibrating data at higher metallicities than shown in Fig- 
ure[6]is far more desirable from the standpoint of using SPS to 
model galaxies since galaxies are generally of higher metallic- 
ity. Unfortunately, old metal-rich star clusters are rare. Clas- 
sic metal-rich clusters include NGC 6791, NGC 188, NGC 
6553, and NGC 6528; the latter two belong to the bulge and 
thus are along sightlines of high extinction. 

In Figure |7] we compare multi-color CMDs of NGC 6791 
and NGC 188 to our SPS predictions. The consideration of 
multi-color CMDs provides a powerful constraint on SPS 
models in part because the reddening value cannot be arbi- 
trarily varied to fit a color without potentially worsening the 
fit at another color. NGC 188 and 6791 are open clusters in 
crowded fields and therefore CMD construction is challenging 
without radial velocities and/or proper motions to isolate true 
members. This complication should be kept in mind when 
considering model comparisons below. 

NGC 6791 is assumed to have an age of 8.5 Gyr, E(B- 
V) = 0.17, [Fe/H1= +0.4, \afFe]= 0.0 , and m-M = 13. 
dCarnev etalj 12001 lOrigliaet all 120061: [Kaliraiet al.l l2007). 
JK and BVI photometry are provided by lCarney et al.1 (2005) 
and IStetson et al . (2003), respectively. The cluster parame- 
ters for NGC 188 are taken to be E(B-V) = 0.09, [Fe/H]= 
-0.06, fa/Fe1= 0.0, m-M = 11.1 and an age of 7 Gyr 
dSaraiedini e t alJ[T999t lCarretta et al.|2009l). UBVRI P hotom- 
etry of this cluster is provided by ISaraiedini et al. (1999); 
near-IR photometry is extracted directly from the 2MASS 
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FIG. 10.- — SSP UV color s as a function of met allicity, for old populations ft ks 13 Gyr). Model predictions are compared to observations of star clusters in M3 1 
(Rev et aljr2007l) . the MW JDorman et"al]| 19951) , and M87 Jgohn et al. 2006). A variety of models are considered to explore the effects of post-AGB (PAGB) 
stars and blue HB (BHB) stars. These models, which are described in detail in the text, are constructed with the Padova isochrones. Model predictions from 
M05, BC03 are al so included, as is a model with the BaSTI isochrones. In addition, the approximate location of Elliptical galaxies from the SAURON survey 
IJeong e t al. 2009) are included in the right set of panels a s erro r bars at log(Z/ZQ) = 0.0. Approximate locations of two metal-rich MW star clusters NGC 188 
and NGC 6791 are shown as open squares (see Sohn et al. 2006). Notice the different y— axes in the left and right panels. 



point source catalog . 

We use the BaSTI isochrones for the following com- 
parisons because those isochrones have been computed for 
log(Z/Z0) = +0.3, which is comparable to the metallicity of 
NGC 6791 and is higher than the highest metallicities avail- 

10 The 2MASS point source catalog was queried us- 
ing the NASA/IPAC Infrared Science Archive, available at 

http : / /ir sa . ipac .caltech.edu/ applications / Gator/ . 



able from the Padova isochrones. Our model CMDs are 
compared to NGC 6791 and NGC 188 in Figure |7j for the 
same reddening values, ages, and distance moduli as quoted 
above, but for slightly difference metallicities. Metallicities 
of log(Z/Z ) = +0.3 for NGC 6791 and log(Z/Z Q ) = +0.0 for 
NGC 188 are adopted based on the available isochrones. 

The agreement between data and models is generally good, 
including the location of the main sequence turn-off, the sub- 
giant branch and RGB, and the position of the red HB. 
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The model fails to reproduce the low mass end of the main 
sequence in NGC 6791 in the sense that the model isochrone 
is on the blue side of the data in B-V and V-K and on the red 
side in V -I. Of course, this disagreement is of little conse- 
quence for integrated properties since low mass stars are faint. 
The disagreement at low masses is probably due to uncertain- 
ties associated with color- T e ff relations for low-m ass stars at 
high metallicity dCassisi et alJl2000tlCass isi 2007). It appears 
that the model RGB V-K color becomes progressively bluer 
than the RGB of NGC 6791 at higher luminosities, although 
it is difficult to make firm statements owing to the paucity of 
data at the bright end. Notice that this discrepancy is in the 
same sense as was seen for integ rated V -K colors of metal- 
poor MW and M31 clusters (see ^4.3. 11 1. 

The model CMDs perform well when compared to NGC 
188. The only obvious discrepancy is in the U-B color of 
the RGB. Thankfully, this is of little consequence for inte- 
grated colors because the integrated U-B color is determined 
primarily by the main sequence turn-off point and sub-giant 
branch, not the RGB (see below). 

For the purpose of modeling entire galaxies, we are inter- 
ested in the ability of SPS models to accurately predict the 
integrated colors of star clusters. It is therefore worth consid- 
ering the extent to which constraints from resolved CMDs can 
constrain integrated colors of clusters. 

In Figure[8]the integrated colors corresponding to the CMD 
of NGC 188 are explored. Integrated colors are computed by 
including all stars fainter than a given apparent V-band mag- 
nitude. It is clear from this figure that integrated blue colors 
such as U-B and B-V are well constrained even at relatively 
faint magnitudes in V. This occurs because integrated blue 
colors are determined primarily by the main sequence turn- 
off point and sub-giant branch. In contrast, redder and near- 
IR colors are highly sensitive to the brightest one magnitude 
(1 1 < V < 12) for the simple reason that the RGB and AGB 
dominate the near-IR light. The brightest one magnitude in 
V contributes ss 20-40% of the integrated V-band light (de- 
pending on whether or not one includes the red clump HB), 
and the V -I and V-K colors redden dramatically in this 
brightest magnitude. Thus, unless the brightest, and rarest 
stars in clusters can be reliably identified and used for model 
calibration, model constraints in CMD-space will not guaran- 
tee accurate model integrated colors. 

4.3.3. Spectral indices 

We now turn to an analysis of the spectral indices of old 
populations over a range of metallicities. Observed indices 
are measured fro m the catalog of optic al spectra for 40 MW 
star clusters from Sc hiavon et alj d200 5). The indices are mea- 
sured in the exact same way for both the observed spectra and 
the models, ensuring reliable comparison. The FSPS model 
spectra are at slightly higher resolution than the MW cluster 
spectra, so the model spectra have been broadened to match 
the resolution of the data. The spectral resolution of the BC03 
model is comparable to the observed spectra, and so the BC03 
model spectra were not broadened. 

Figure|9]presents the spectral indices D„4000, HSa, Fe5270, 
and Mgb as a function of metallicity for the Schiavon et 
al. sample and the FSPS and BC03 SPS models. The 
FSPS predictions utilize the Miles empirical stellar library, 
while BC03 adopt the STELIB empirical library. For these 
ages the Padova predictions are very similar to the BaSTI 
isochrones except for HSa at low metallicity because the 
BaSTI isochrones have a hotter HB. This figure also demon- 



strates the importance of BS stars and an extended/bl ue H B 
by including models with 5bs = 2 and /bhb = 0.5 (see t j3.ll ). 

There are several important conclusions to draw from this 
figure. First, the models are able to reproduce the observed 
metallicity dependence of the Mgb and Fe5270 indices over 
the full observed range. At low metallicities this success is 
driven largely by the fact that the bulk of the metal-poor 
stars in the empirica l stellar library contains a-enhanced stars 
dMilone et alj|2009t) . While the isochrones used in FSPS are 
solar-scaled 11 , recent work has shown that for old popula- 
tions the effects of a-enhancement are strong on the stellar li- 
braries and weak on the isochrones ( Schiavon 20071 iLee et alj 
2009cE). 

The default SPS models (i.e., without BS or blue HB stars) 
predict D„4000 in excess of observed values at solar metallic- 
ity, and the same models predict HSa strengths too weak com- 
pared to the observations at both very low and solar metallic- 
ity. The disagreement at matching HSa at very low metallicity 
is a consequence of the well-known fact that most low metal- 
licity clusters have a blue HB. Indeed, the models do a better 
job of matching the observed HSa strengths at low metallicity 
for /bhb > 0.5. 

The failure of the models for both D„4000 and HSa at so- 
lar metallicity is more striking, and has potentially signifi- 
cant consequences for understanding the spectra of metal-rich 
galaxies. The two solar metallicity clusters are NGC 6553 and 
NGC 6528; both are in the bulg e. The CMD of NG C 6553 
shows an abundance of BS stars dZoccali et alJl200Tb . and so 
the addition of such stars to the models in Figure[9]is justified. 
Indeed, adding BS stars results in much better agreement in 
HSa and moderately better agreement for D„4000. The addi- 
tion of blue HB stars also a lleviates much of the discrepancy 
between models and data. ILee et al] d2009bl) compare their 
model to the same dataset and find that for solar-scaled mod- 
els the HSa strengths are too weak compared to the observed 
Z ss Z© clusters, echoing our results regarding our default SPS 
model. 

Main sequence turn-off stars are the dominant stellar phase 
controlling the strength of D„4000. The effective temperature 
of turn-off stars would have to be increased by ~ 1000K in 
order to alleviate the discrepancy at solar metallicity. Such an 
increase would yield a CMD in gross conflict with observa- 
tions. Modifications to the standard stellar phases thus cannot 
reconcile the data and model D„4000 strengths. Either BS, 
blue HB stars or both are therefore the most likely candidates 
for explaining the discrepancy. 

4.3.4. Ultraviolet colors 

In the absence of evolved or young stars, the ultraviolet is 
sensitive to the temperature of main sequence turn-off stars. 
Since the temperature of turn-off stars depends on metallicity 
at fixed age, ultraviolet colors such as NUV - V and FUV - V 
will, in the absence of evolved or young stars, redden substan- 
tially with increasing metallicity. 

In reality, the interpretation of ultraviolet colors of old pop- 
ulations is complicated by the existence of hot, evolved stars, 
including post-AGB, extended/blue HB, and BS stars (e.g., 
i Burstein et al.ll 1 988fc iGreggio & Renzinilll990t iDorman et alj 
Il995t IO'Connel]|ll999t IBrown et alJ 1200a IHan et alJI2007E 
It is beyond the scope of this paper to consider in detail the 
constraints on the relative importance of these various popu- 

1 ' Although recall that the BaSTI database does provide a- enhanced stel- 
lar evolution models. 
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FIG. 1 1. — Color-color diagrams comparing observed luminous red sequence galaxies to both our FSPS model with the modified Padova isochrones, and the 
models of BC03 and M05. The contours enclose 40%, 68%, and 95% of the observed galaxies. For each model, the colors of dust-free SSPs are plotted for ages 
5 < t < 13 Gyr for three metallicities: \o°(Z/Zq) ks -0.2,0.0,0.2 (corresponding to three lines for each model). Symbols mark the location of t = 13 Gyr for 
each metallicity. In the lower right panel, all model predictions are shown as symbols for clarity. 



lations afforded by ultraviolet colors. Instead, our goal will be 
to demonstrate that plausible variations in the relative weights 
assigned to post-AGB and blue HB stars can explain the ul- 
traviolet properties of observed star clusters. 

Ultraviolet colors of star clusters in M31 dRev et alj|2007l) . 
the MW (as com piled by iDorman et all 11995b . and M87 
(ISohn et al.1 120061) are compared to SPS models in Figure 
[TOl The approximat e location of E lliptical galaxies from the 
SAURON survey dJeong et alJ l2009) are also included in the 
figure. We include model predictions from our FSPS code 
for both the Padova and BaSTI isochrones, along with model 
predictions from BC03 and M05. Recall that M05 does not 
include post-AGB stars and so the FUV photometry, which is 
sensitive to such hot stars, cannot be reliably compared to the 
data. M05 does include blue HB stars for metal-poor models. 
BC03 includes post-AGB stars but does not include blue HB 
stars at any metallicity. 

For the FSPS predictions with the Padova isochrones, we 
include several model variations. In particular, we consider 
models 1) without post-AGB and without blue HB stars; 2) 
with post-AGB but without blue HB stars; 3) without post- 
AGB stars and with a blue HB fraction of 10%; 4) with post- 
AGB stars and a blue HB fraction of 100%; 5) with post-AGB 
stars and a blue HB fraction that varies with metallicity as 
/bhb = -O.751og(Z/Z )-O.5 for log(Z/Z ) < -0.7 and zero 
otherwise. 

Before considering the insights gained from model compar- 
isons, it is worth drawing attention to the considerable vari- 



ation in observed ultraviolet colors of star clusters at fixed 
metallicity. NUV - V colors vary by w 0.5 mags at fixed 
Z, while FUV — V varies by > 1 mag. This phenomenon is 
known as the 'second-param eter' effect, an d has bee known 
for decades (see the review in Catelan 2009). Parameters that 
might be responsible for the additional variation at fixed Z 
include age, helium content, mass-loss, stellar rotation, and 
magnetic fields. The particular solution is not as important as 
the observed variation. It is difficult to imagine constraining 
an SPS model in the UV given these variations, and it will 
therefore continue to be very challenging to extract useful in- 
formation from the UV flux of quiescent galaxies. 

The locus of ultraviolet colors of Elliptical galaxies do not 
correspond to star clusters of any metallicity in either M31, 
the MW, or M87. The NUV -V colors correspond to metal- 
rich clusters, while the FUV - NUV colors are similar to 
metal-poor clusters, and the FUV - V colors of the Ellipti- 
cal sample are significantly redder than any star cluster. This 
conflict may be due in part to the fact that galaxies are not 
mono-metallic. The ultraviolet is quite sensitive to metallic- 
ity, and so accurate modeling of the metallicity distribution 
function is essential for these objects. These and related is- 
sues regarding s tellar pop ulation modeling of Ellipticals will 
be discussed in g3S21 

It should also be noted that the star clusters in M87 are sig- 
nificantly bluer in ultraviolet colors at fixed metallicity com- 
pared to clusters from the MW and M3 1 . The source of this 
difference is not known, although differences in cluster age 
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FIG. 12. — Comparison of spectral indices for a sample of SDSS galaxies (68% and 95% contours) to our FSPS model predictions (lines). Here, our FSPS 
model uses the empirical Miles stellar library and the Padova isochrones. SDSS galaxies are divided by a > 200km s (solid contours) and 100 < a < 200km s 
(dashed contours). Model predictions are displayed for dust-free SSPs with metallicities of Zq (solid lines) and 1.5Zq (dashed lines) for 2 < t < 13 Gyr. Model 
spectra are Doppler broadened to cr = 200kms~' . In the upper left panel, symbols denote ages of 2.2.3.5,5.6, 8.9, 14.1 Gyr, in order of increasing Z)„4000. 



or helium abundance are possible explanations (So hn et all 
I2006t iKavirai et al.ll2007l) . Nevertheless, it is clear that even 
if models can be constructed that match the MW cluster data, 
such models may not be applicable to other systems. 

Several conclusions can be drawn from the comparisons 
between models and data. First, a 10% blue HB fraction at 
all metallicities produces ultraviolet colors that are compa- 
rable to post-AGB stars, demonstrating the degeneracy be- 
tween these two evolutionary phases when considering UV 
colors. The BaSTI isochrones and the Padova isochrones 
with variable blue HB fraction are both able to fit a major- 
ity of the metallicity-dependent ultraviolet data. At metal- 
licities log(Z/Z Q ) > -0.7 these models fail to reproduce the 
observed FUV fluxes, whereas models with 100% blue HB 
fractions can reproduce this data. 

When modeling the ultraviolet flux of galaxies, a flexible 
treatment of hot evolved stars is clearly warranted given our 
inadequate theoretical understanding of the observed ultravi- 
olet colors of clusters, and given the considerable scatter in 
observed ultraviolet data at fixed metallicity. 

4.4. Constraints from quiescent galaxies 

We now consider constraints on SPS models provided by 
massive red sequence, (i.e., quiescent), galaxies. Such galax- 
ies are thought to be composed primarily of old (f > 10 Gyr), 
metal-rich stars. 



4.4.1. Colors 

Figure QT| compares galaxies to models in color-color 
space. We include the FSPS model predictions with the modi- 
fied Padova isochrones in addition to the BC03 and M05 mod- 
els. 

All models predict similar ugr colors. One clearly sees 
the almost perfect degeneracy between age and metallicity in 
this color-color space. It is also clear that all models predict 
ugr colors that are substantially redder than the data, assum- 
ing th at such galaxies are composed of old, metal-rich stars 
(e.g., Iwbrthev et al.l 11992k iTrager et al.1 12000t iThomas et al.l 
120051: iGraves & Schiavon 2008^ The discrepancy between 
observed and predicted g—r colors is > 0.1 mag s. Potential 
remedies to this discrepancy are explored in H4.5\ 

The models perform reasonbly well in the near-IR, al- 
though there is a notable disagreement in the J — H color for 
all models. In the r-z vs. H-K plane the FSPS and M05 
models perform well, while BC03 fits the data only for solar 
or sub-solar metallicities. The BC03 and M05 models pre- 
dict strong variations in Y - K with metallicity, in contrast to 
FSPS. Given the narrow distribution of observed Y—K col- 
ors, the BC03 models would predict nearly all red sequence 
galaxies to be solar metallicity, while M05 would predict that 
all were super-solar. In contrast, FSPS allows for a range 
of metallicities from sub- to super-solar, in better agreement 
with detailed analysis of the spectra of red sequence galaxies 
(e.g.. IThomas et al.ll2005l) . 
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FIG. 13. — Spectral indices for a sample of SDSS galaxies with a > 200kms~' (contours), compared to several models for dust-free SSPs. FSPS model 
predictions are shown for both the Padova and BaSTI isochrones where in each case the empirical Miles stellar library is used in the model construction. 
Predictions from the BC03 model are also included. Models are shown for metallicities of Zq (solid lines) and 1 .5Zq (dashed lines) and for ages in the range 
2 <t < 13 Gyr. Model spectra are Doppler broadened to a — 200kms~' . 



Perhaps the most important conclusion to be reached from 
these color-color comparisons is that no model is able to re- 
produce the observed locus of red galaxies for all near-IR col- 
ors, and in addition no model can match the locus of galaxies 
in the u—g vs. g—r plane. 

It should also be noted that within the context of one model, 
optical and near-IR color combinations, in particular g — r vs. 
Y—K, are in principle capable of breaking the age-metallicity 
degeneracy. In practice it is not currently possible to use such 
colors to break this degeneracy for the twofold reason that the 
models cannot predict the observed locus of galaxies in that 
space an d the models do no t compare well wi th each other 
(see also lEminian et alj|2008t ICarter et al.ll2009l) . 

The red and near-IR colors of these galaxi es are determined 
primarily by the RGB. As noted in 34.3.21 resolved CMDs 
do not help constrain model predictions for integrated red 
and near-IR colors because these colors are determined by 
the brightest and rarest RGB stars. Homogeneous, integrated 
near-IR photometry for a sample of old metal-rich clusters 
would be invaluable for constraining models in this regime. 

4.4.2. Spectral indices 

Figure[T2]compares spectral indices of SDSS galaxies to the 
FSPS model constructed with the Padova isochrones and the 
Miles empirical stellar library. The data are split according 
to their velocity dispersions. All models in this section are 
broadened to a = 200 km s" 1 in order to compare to the more 



massive data bin. Model predictions are shown for dust-free 
SSPs at Z = Z Q and Z = 1 .5Z Q for ages 2 < t < 13 Gyr. 

The index D„4000, which measures the strength of the Call 
lines and the balmer decrement, is often taken to be sensi- 
tive primarily to age, and this is borne out in the figure. The 
majority of the other indices are sensitive to a combination 
of age and metallicity. The index C4668, which is a blend 
of many metallic line s including Fe, Ti, Cr, Mg, Ni, and C2 
(Wort hey et al.l [l994). is the most sensitive to metallicity of 
the indices shown, at fixed D„4 000. The index [MgFe]' was 
defined by iThomas et al.l d2003l) in order to be insensitive to 
a-enha ncement. The G4300 index is sensitive to both CH 
and Fel (Worth ey et al.llT994b . H<5>t is often used to constrain 
the fraction of r ecent star formation bec ause of its sensitivity 
to A stars (e.g.. iKauffmann et al.ll2003h . and so the reliabil- 
ity of SPS models in the HSa - D„4000 plane is particularly 
important. For old populations, the interpretation of H5a is 
complicated by the presence of nearby absorption due to the 
CN molecule. This implies that the HSa index can be sensi- 
tive to non-solar abund ance ratios, depe nding on the precise 
definition of the index dProchaska et al.l 120071) . The trend is 
such that increasing a-enhancement will result in a weaker 
measured HSa index because the flux redward of the line will 
be depressed due to increased CN absorption. 

Figure [13] compares several SPS models to the a > 
200kms _I SDSS sample. The BC03 model predictions are 
compared to the FSPS models using both the modified Padova 
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FIG. 14. — Comparison of the spectral indices D„4000 and HSa between 
SDSS galaxies with a > 200kms~' (contours) and various models con- 
structed with FSPS adopting the Padova isochrones (as in previous figures). 
Model spectra are Doppler broadened to a = 200 km s^ 1 and span the age 
range 2 < t < 13 Gyr, with the oldest point marked by a solid symbol. Blue 
HB (BHB) stars, young stars (with an age of 1 Gyr), metal-poor stars, and 
BS stars (with Sgs = 2) are considered as additions on top of a base dust-free 
SSP of either solar or super-solar metallicity. The locations of two ~ Zq 
MW bulge star cluste rs are included (crosses) from the integrated spectra of 
ISchiavo n et al. (2005). The MW star cluster spectra were Doppler broadened 
to a — 200kms~' for comparison to the galaxy data. 



and BaSTI isochrones with the Miles empirical stellar library. 
The M05 model is not included because the spectral library in 
that model (The Kurucz library) is not of sufficient resolution 
to measure these indices. 

The BaSTI and Padova isochrones have very similar main 
sequence turnoff points at old ages (see Figure [TJ. The pri- 
mary difference between these isochrones is that the BaSTI 
RGB is approximately 100K hotter than the Padova RGB. 
Blueward of A = 5000A the main sequence and sub-giant 
branch dominate the flux, and so one would expect very little 
differences in spectral indices between these two isochrone 
sets in the blue, for a fixed set of empirical stellar spectra. 
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FIG. 15. — Comparison of the observed ugr colors of luminous red se- 
quence galaxies (contours) to FSPS models. All models are dust-free SSPs 
at Z = Zq for 5 < t < 1 3 Gyr, with symbols marking the location of t = 1 3 
Gyr. The standard model is compared to models where metal-poor stars, BS 
stars (with Sbs = 2), and young stars with an age of 1.0 Gyr, are added on top 
of the standard model in the proportions noted in the legend. 



Indeed, in Figure [13] there are almost no differences between 
the two isochrone sets for those indices that are measured at 
A < 5000A. The only index that shows appreciable differ- 
ences is Ti02, which is measured at A « 6230A, where the 
RGB contributes more to the integrated flux than the main se- 
quence turnoff point. 

The BC03 and FSPS models compare favorably for the 
D„4000, [MgFe]', G4300, and Ti0 2 indices (except for the 
Z = 1 .5Z Q models for the latter index). The Z = 1 .5Z Q FSPS 
model predictions for C4668 compare more favorably with 
the data than the BC03 model of comparable metallicity, 
while the converse is true for the Fe4531 index. The dif- 
ferences between the BC03 and FSPS models are most dra- 
matic in the H6 A -D„4000 plane, where the BC03 HS A EWs 
are larger than the FSPS predictions by « 1.3A, and are in 
significant discord with the data. Thi s discrepancy has been 
noted previously bv lWild et al.1 d2007l) . and is most likely due 
to deficiencies in the STELIB empirical stellar library used by 
BC03. 

One of the most striking results of these comparisons is 
that all models tend to over-predict the D„4000 strength and 
under-predict the HSa strength, especially for Z = 1 .5Z©. As 
can be gleaned from the symbols in the upper left panel of 
Figure [12] at Z = 1 .5Z the models reproduce the observed 
D„4000 locus of massive galaxies for ages 3-6 Gyr, while for 
Z = Zq the observations are reproduced for somewhat older 
ages of 6-10 Gyr. F rom a detailed analy sis utilizing mul- 
tiple spectral indices, Thom as et al.l (|2005) find most mas- 
sive ellipticals to have solar to super-solar metallicities and 
ages > 10 Gyr (see also, e.g.JBower et al.lll992tlBender et alJ 
U996tlZiegleret al]ll999tlGraves & Schiavonll2008l) . In light 
of these results, an explanation is required for the discrepancy 
between observed and predicted D„4000 and HSa values. Pos- 
sible explanations are discussed in the following section. 

Before considering possible explanations we point out that 
the proper inclusion of a-enhancement will only exacerbate 
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the discrepancy between models and data in regards to the 
HSa index. This occurs bec ause q-enhancement w ill result in 
an even weaker HSa index (Proch aska et al.ll2007l) . 

4.5. Resolving the model-data discrepancies for quiescent 

galaxies 

In the previous section two striking disagreements were 
noted between SPS models and a sample of massive red se- 
quence galaxies. Assuming that such galaxies are composed 
of 'normal' stars 12 of at least solar metallicity and with ages 
> 10 Gyr, all SPS models predict ugr colors that are too red, 
D„4000 strengths that are too large, and HSa strengths that are 
too weak. 

There are two possible classes of solutions to this discrep- 
ancy. In the first class, an appeal is made to deficiencies in 
either the stellar atmospheres or stellar evolution calculations. 
Within this class, an appeal to stellar evolution means con- 
sidering inadequacies in standard phases such as the main se- 
quence and its turn-off point, and the RGB. Solutions of the 
second class abandon one or more of the assumptions regard- 
ing the composition of such galaxies. In this second class ap- 
peals are made to minority populations of either young stars 
or metal-poor stars, or exotic aspects of stellar evolution in- 
cluding blue stragglers and mass-loss leading to an extended 
horizontal branch. These classes are conceptually different in 
that the former appeals to a modification of the bulk metal- 
rich old population while the latter appeals to additions of mi- 
nority populations. We consider both of these classes below. 

4.5. 1. Problems with stellar atmospheres and stellar evolution 

We find that the discrepancies are unlikely to be due to the 
stellar spectral libraries. We have considered both the theo- 
retical BaSeL library and the empirical Miles library and find 
similar disagreements in the colors. The empirical STELIB 
library used in BC03 yields similar disagreements between 
both the g—r color and the D„4000 strengths, further suggest- 
ing that the libraries are not at fault 

These results are at odds with Marast on et al.l (120091) who 
showed that the optical spectra and g — r colors of SSPs at Z Q 
are markedly different when comparing models cons tructed 
using the BaSeL spectral library and the empirical iPicklesi 
(1998) library. The purported discrepancies caused by the 
spectral libraries cannot be reproduced with our FSPS code 
when using the empirical Miles library. To test that this was 
not due to differences between the Miles and Pickles library, 
we also compared spectra generated from the BC03 model 
using both the BaSeL and Pickles library and again find no 
substantial difference in g — r colors. We have also compared 
the spectra of individual stars from the Pickles, Miles, and 
BaSeL libraries for the same stellar parameters discussed by 
Maraston et al. and again cannot reproduce their results (i.e., 
their Figure 4). 

We also point out that in the context of the explanation pro- 
posed by Maraston et al., the BaSeL library should provide a 
better fit to the data when the colors are k-corrected to z = 0.0 
(as opposed to z = 0.1), because their discrepancy arises in a 
region of the spectrum sampled by the r-band redshifted to 
z = 0.1. However, when we compare the colors of z < 0.03 
red sequence galaxies we find that the model discrepancies 
remain. 

12 'Normal' here is meant to exclude such stellar phases as blue stragglers 
and hot horizontal branch stars. 



The discrepancy is also unlikely to be due to unresolved 
issues with canonical phases of stellar evolution such as the 
main sequence and RGB, because the aspects of stellar evo- 
lution relevant for ugr colors are adequately described by the 
models for old, metal-rich star clusters. This claim is based 
on two observations. First, from Figure|7]it is clear that FSPS 
can reproduce the main sequence turn-off point and base of 
the RGB of old, metal-rich star clusters. This is especially the 
case for UBV colors of the clusters NGC 6791 and NGC 188. 
To the extent that there are residual discrepancies, the models 
tend to be somewhat too blue compared to the data, which is 
in the opposite sense of the discrepancies noted for galaxies. 
Second, in Figure[8]it was shown that for blue colors such as 
B-V, one need not probe the brightest portions of the RGB 
in order to achieve convergence in the integrated color (the re- 
sults for g — r are very similar to B-V). Therefore, integrated 
g—r colors of old metal-rich clusters should be accurate be- 
cause the main sequence turn-off and the base of the RGB are 
relatively well-fit by the model in the CMD. 

In sum, we believe that the discrepancies between model 
and data ugr colors and spectral indices are not a reflection 
of the limitations of either stellar atmospheres or canonical 
aspects of stellar evolution. 

4.5.2. Minority populations in predominantly old metal— rich 
galaxies 

There are two remaining possible solutions to the disagree- 
ment between models and data for massive red sequence 
galaxies, which we lump into one class because they share 
the common property of being minority additions to the dom- 
inant, old metal-rich population. The first solution consists of 
modifications to the SSPs. This solution includes the addition 
of BS or blue HB stars. 

The second solution to these disagreements lies not in try- 
ing to fix the SSPs but rather in relaxing the implicit assump- 
tion that these massive red galaxies are composed of mono- 
metallic, coeval stars. Indeed, estimates of the metallicities of 
stars in our own and nearby galaxies, either from direct spec- 
troscopy or via the colors of giant branch stars, always yield a 
metallicity distribution function tha t is approximately lognor- 
mal with a full width of ~ 1 dex jRichll 1 990c iGrillmair et alJ 
1996b lHarris et al.l 11999b IZoccali et all 12003b IWorfhev et alJ 
20051 ISarajedini & Jablonkal2005blKoch et alj2006l) . Models 
for the chemical evolution of galaxies also generically pre- 
dict that galaxies should contain stars of a range of metal- 
licities, including a non-trivial frac tion of metal-poor stars 
(ISearle & SargentfT9 72: Pag el 19971) . 

Extreme/blue HBs exist in the solar and super-s olar metal- 
licity s tar clusters NGC 188 dLandsman et alll998h and NGC 
6791 dLandsman et all 11998b IKalirai et alJ [20071). in addi- 
tion to moderately sub-solar clusters dRich et aLlll997h . Re- 
solved UV photometry in th e nucleus of M32 ha s shown that 
~ 5% of HB stars are hot dBrown et all 120001) . Moreover, 
man y metal-rich clusters contain copious BS populations 
(e^g jLandsman et alJ[T998l: IZoccali etalj|200lt IKalirai et al J 
2007). These observations suggest that it would be prudent 
to consider models that allow for blue/extended HBs and BS 
stars to be associated with metal-rich populations. In our 
treatment we are agnostic as to which population (metal-rich 
or metal-poor) these hot stars are associated with. 

These various effects are explored in Figures [14] and Q3] In 
these figures we explore the impact of various modifications 
to the base dust-free SSP models. Modifications include: 1) 
a small fraction (1-3%) of young (f = 1.0 Gyr) stars; 2) a 
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FIG. 16. — Color-color diagrams comparing observed K+A (post-starburst) galaxies (grey symbols) to both our FSPS model and the models of BC03 and M05. 
For BC03 and M05, the colors of SSPs are plotted for ages 0.8 < t < 2.5 Gyr for three metallicities: log(Z/Zo) ~ —0.2,0.0,0.2 (corresponding to three lines 
for each model). FSPS models are shown for a r— model SFH with r = 0. 1 Gyr, for the same age and metallicity ranges as the other models. The choice of a 
composite SFH for FSPS was chosen for cosmetic reasons only. Symbols a long the model predi ctions mark t = 2.5 Gyr. In the left panel the arrow indicates 
the change in colors that would occur from adding dust with Ay = 0.5 and a Calzetti et al. 1 1994) reddening curve. As discussed in the text, the addition of an 
underlying old population moves the models along the reddening vector, which is also along the locus of observed griz colors. 



population of blue HB stars (/bhb = 0.2); 3) a population of 
BS stars (Sbs = 2); 4) a minority population (5%) of metal- 
poor stars. 

All of these modifications move the base SSP models in the 
correct qualitative sense; that is they result in lower values of 
D„4000 and larger HSa strengths, and bluer ugr colors. These 
trends result from the fact that each modification is adding 
a population of relatively hot stars, whether in the form of a 
hotter main sequence turnoff (as happens when metal-poor 
stars are added), hot evolved stars including BS and blue HB 
stars, or hot young stars. Hotter stars have lower values of 
D„4000, larger balmer absorption, and bluer continua. 

4.5.3. Summary and outlook 

The discrepancies discussed above between SSP predic- 
tions and massive Elliptical galaxies h ave been known, in 
one form or another, for dec ades (e . g.. iFaber 119771; lOconnelll 
[Q iRosd 11985b ISchweizer et al l 
1992t IWorthev et al.1 119961: 



1980; Gunnetal. 1981 
1 990t ISchweizer & Seitzerl 

Marast on & Tho mas 200oT lCaldwell et al.ll20d3T lTrager et al. 
20051: ISchiavonl2 007 ) . Our contribution to this issue has been 
twofold: 1) to highlight the connection between the indices 
and broadband ugr colors, in the sense that whatever hot star 
population is adopted, it should alleviate the discrepancies 
seen in both indices and colors; 2) to remind the user of 
SPS models that reliable results concerning the ages and 
metallicities of massive Ellipticals cannot be attained unless 
these minority hot star populations are carefully considered. 

Many previous studies have considered one or more of 
the minority populations considered here in as a remedy to 
the model discrepancies. Schiavon (2007) disfavors appeals 
to blue HB stars based on the consideration of higher-order 
balmer lines in Ellipticals. However, this conclusion rests on 
the assumption that the temperature distribution along the HB 
in other galaxies is similar to that seen in the MW. While this 
may seem to be a relatively benign assumption, recent UV 



observations of clusters in M87 cast doubt on this ( So hn et al.l 
120061) . If the M87 result holds up to further scrutiny, it will 
mean that insights from the MW star cluster population can- 
not be readily applied to other galaxies. This would be a 
dramatic setback to our attempts at understanding the stellar 
populations in other galaxies, and would have potentially pro- 
found implications for our understanding of stellar evolution. 

The addition of metal-poor stars is rarely considered with- 
out also considering the addition of blue HB stars. In 
this work we have separated these two populations, since 
the observed relation between metallicity and HB morphol- 
ogy in the MW is com plex and displays considerable scat- 
ter dPiotto et alJ 120021) . It should be remembered that in 
every galaxy where metallicity distribution functions have 
been measured, there exist metal-poor s tars (log(Z/Zp) <-, 
— 1) in sufficient numbers (~ 1-5%; e.g., [Harris et al l ll999t 
IZoccali eTa l. 2003; Sarai edini & Ja blonka 2005) to effect the 
integrated blue/UV flux. Given this observational fact, metal- 
poor stars must be included in the models; the only question 
that remains is in what proportion. 

It is something of a disappointment that the identification 
of the hot star population in Ellipticals remains unresolved 
more than three decades after it was first discussed. Despite 
enormous advances in both the SPS models and the quality of 
the data, the fundamental problem remains that these various 
hot star populations occupy essentially the same (or similar) 
ranges in r e ff. We briefly comment on several possible future 
directions that might help separate the relative influence of 
these various populations. 

If young stars or metal-poor stars are the dominant contrib- 
utor, then because their metallicities will differ from the bulk 
population, metallicities measured in the blue/UV will dif- 
fer from metallicities measured in the red. Models will have 
to be constructed to determine whether or not this effect is 
measurable. If it is then this could be a promising means of 
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discriminating between the various hot star candidates. 

If young stars are the culprit, then lookback studies of El- 
lipticals will uncover Ellipticals with younger and younger 
minority populations. Such a study would suffer from the 
typical difficulties in lookback studies, namely that pre- 
cise descendant-progenitor relations are difficult to construct. 
However, with new and upcoming wide-field spectroscopic 
surveys we believe that the difficulties associated with look- 
back studies can be mitigated. Lookback studies might also 
be ab le isolate the importance of blue HB stars ( Atl ee et al.l 
l2009h . 

A final, somewhat more speculative possibility, is the use of 
SBFs in the blue to discriminate between these various pop- 
ulations. While these populations occupy approximately the 
same range of r e ff, they all occupy rather different ranges in 
Lboi, and they will thus produce different SBF magnitudes in 
the blue. We are currently constructing models aimed at de- 
termining whether or not such an effect is observable with 
current facilities. 

4.6. Constraints from post-starburst galaxies 

K+A, or post-starburst, galaxies offer a unique constraint 
on intermediate age populations. The prominent A star spec- 
tra in these galaxies without any accompanying Ha emission 
(i.e., no current star formation), indicates that they are domi- 
nated by stars with ages wl-2 Gyr. Such galaxies can there- 
fore provide potentially strong constraints on the importance 
of TP-AGB stars, since such stars peak in prominence at ages 
w 1-2 Gyr. 

Fi gure [TBI compare s the colors of K+A galaxies defined by 
Ouint ero et al.l (120041) to the colors of dust-free SSPs at ages 
0.8 < t < 2.5 Gyr for metallicities log(Z/Z Q ) « -0.2, 0.0, 0.2. 
Model predictions from FSPS (using the Padova isochrones), 
M05, and BC03 are shown. The FSPS model predictions are 
shown for a r-model SFH (SFRoc e~'l T ) with r = 0.1 Gyr 
for cosmetic reasons only. Recall that the TP-AGB phase 
in the Padova isochrones was modified in order to produce 
acceptable agreement with both the MC star clusters and these 
post-starburst galaxies (although the modificat ions a re slight 
for intermediate ages at solar metallicity; see £|3 . 1 . 3b . While 
not shown, FSPS predictions using the BaSTI isochrones are 
similar to the Padova isochrones. A reddening vector is shown 
to indicate the shift in colors expected f or a uniform sc r een of 
dust with optical depth A v = 0.5 and a ICalzetti et al.l (119941) 
attenuation law. 

The ratio of A to K stars in the observed spectra decreases 
with increasing g — r and i - z colors, indicating that the se- 
quence in the left panel is primarily a sequence in age, and 
not dust nor metallicity. This observational fact will be a 
significant discriminator amongst the SPS models considered 
herein. 

The FSPS model predictions encompass the range of ob- 
served i—z, g — r, and Y — K colors. The dominant variable 
driving the change in model colors is age, in agreement with 
the data. The BC03 model predictions are somewhat too blue 
in i - z, although we believe that the agreement is acceptable 
given the fact that post-starburst galaxies are composite pop- 
ulations whereas we have modeled them as single populations 
(see below). 

In contrast to the FSPS and BC03 models, the M05 model 
performs poorly in g— r and i—z colors. In particular, model 
sequences of increasing age produce trends in these colors 
that are orthogonal to the observed trends with age. In the 
context of the M05 model the increasing g—r and i—z colors 



would be interpreted as a sequence of either increasing dust 
content or increasing metallicity. This interpretation would be 
in contrast with the observed variation in the spectral proper- 
ties. It is also difficult to explain the bluer half of the observed 
K+A galaxies in griz colors unless unrealistically low metal- 
licities are considered. Appealing to younger ages does not 
remedy this either, as M05 model predictions at younger ages 
(t < 0.8 Gyr) continues to diverge from the locus of observed 
colors. 

As discussed in previous sections, it is difficult to reliably 
calibrate SPS models against galaxies because galaxies con- 
tain stars of a range of metallicities and ages, and the effects 
of dust cannot be ignored. Here we have been considering 
a dust-free, mono-metallic, coeval population model. For 
the comparisons in this section, these complications turn out 
to have minimal consequences regarding the model compar- 
isons. For example, as shown in Figure [16] the dust redden- 
ing vector lies parallel to the locus of data, and so dust effects 
cannot move the models onto or off of the observed locus, but 
only along the locus. Since our conclusions do not rely on the 
detailed location of the models along the observed locus, our 
results are insensitive to dust effects. 

We have also considered the effects of adding an old pop- 
ulation on top of the intermediate age, post-starburst popu- 
lation. Even in the extreme case where 90% of the mass is 
contained in an old population and only 10% spans the in- 
termediate age range of 0.8 < t < 2.5 Gyr, the red end of 
the model predictions only shift slightly toward redder col- 
ors. The primary effect of an addition of an old population is 
a narrowing of the range in color of the model predictions. In 
any event, the addition of old stars has an effect similar to dust 
attenuation — namely, the models shift along the observed lo- 
cus. We conclude that, in the case of post-starburst galaxies 
and for the colors we have considered, composite population 
effects do not impact our conclusions. 

5. MODEL EVALUATION 

In the previous section we explored the constraints on sev- 
eral SPS models provided by star cluster data in the Local 
Group and M87 and by post-starburst and massive red se- 
quence galaxies. We now summarize these results and pro- 
vide a general evaluation of the models. 

Star clusters in the MCs provide valuable constraints on 
models at intermediate age and sub-solar metallicities. The 
M05 model performs poorly when compared to the extant 
data, while the BC03 model performs better, although the 
BC03 model is somewhat too blue in the near-IR. In M05, the 
colors are too red and the age-dependence is incorrect. The 
incorrect age-dependence in M05 cannot be attributed to the 
details of MC cluster age estimation. Our own FSPS model 
with Padova isochrones is able to achieve excellent agree- 
ment with the data only after we modify the TP-AGB phase 
of the input Padova isochrones. The BaSTI isochrones in 
FSPS fare well, although less-so in the near-IR, which plausi- 
bly reflects their simplistic TP-AGB treatment, and, at young 
ages, which may be due to our use of the BaSTI isochrones 
that suppresses convective overshooting. Similar conclusions 
are reached regarding the FSPS model when considering con- 
straints afforded by SBFs of MC star clusters. 

The MW and M31 provide a large sample of old, metal- 
poor star clusters, but only a small handful of old metal- 
rich clusters. The FSPS (using both the Padova and BaSTI 
isochrones), BC03, and M05 models all make similar predic- 
tions for the UBVRIJHK photometry of old metal-poor sys- 
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terns. The models compare favorably in UBVRI colors once 
uncertainties in the filter transmission curves are considered. 

In the near-IR the models perform less well, where for ex- 
ample V—K colors are too blue by w 0.2 mags at log(Z/Z ) > 
—1.0. The origin of this discrepancy is unclear, but may ei- 
ther be due to an RGB that is too hot (lowering the RGB by 
200K would alleviate this tension) or stellar libraries that are 
inadequate in the near-IR. At higher metallicity, multi-color 
CMDs in UBVIK colors are well reproduced by the FSPS 
model in regards to the location of the turn-off point, sub- 
giant branch, and the lower portion of the RGB. We demon- 
strate that agreement in CMD-space between data and models 
does not imply that the integrated colors are adequately con- 
strained because the integrated colors in the red and near-IR 
are dominated by the brightest stars in the cluster. These stars 
are rarely included in CMD studies because they are rare and 
often saturate the photometry. In any event, even if these stars 
are considered, owing to their rarity it is challenging to reli- 
ably calibrate models in this regime in CMD-space. It is for 
this reason that we focus our calibrations on integrated cluster 
properties, rather than resolved CMDs. 

The spectral indices D„4000 and H5a of MW star clusters 
are not adequately reproduced by either the FSPS or BC03 
model at high and low metallicities, respectively (the M05 
model does not contain sufficient spectral resolution to make 
a robust prediction). These inadequacies are removed if BS or 
hot HB stars are added. 

In the ultraviolet, the FSPS model can reproduce the ob- 
served trends because the relevant phases — post-AGB and 
HB — are handled flexibly and so can be tuned to match the 
data. The BC03 and M05 models generally fare less well, 
either because of the omission of blue HB stars (BC03) or 
post-AGB stars (M05). 

Turning to galaxies, all models predict ugr colors of red se- 
quence galaxies that are too red and the FSPS and BC03 mod- 
els predict D„4000 strengths that are too strong (again, the 
M05 model does not contain the requisite resolution), under 
the assumption that such galaxies are composed solely of old 
stars with Z > Z . These issues are likely closely connected 
to the model failures at reproducing the spectral features of 
metal-rich MW star clusters. These discrepancies can be re- 
moved by adding a minority population of hot stars, coming 
either from metal-poor turn-off stars, hot HB stars, BS stars, 
or young (~ 1 Gyr) stars. 

The FSPS and BC03 models are able to reproduce the op- 
tical and near-IR colors of post-starburst galaxies, while the 
M05 model performs less well. These galaxies likely have a 
preponderance of TP-AGB stars, and so they provide valu- 
able constraints on this phase. 

In summary, our FSPS model is able to reproduce the great- 
est range of observations, due in large part to the flexibility of 
the model. BC03 fares somewhat less well, and M05 performs 
worst of the three. Even where calibrating data is available, 
significant discrepancies remain. Important regions of param- 
eter space do not contain sufficient calibrating data, and so 
the models are even less reliable there. Finally, when model- 
ing predominantly old, metal-rich populations, care must be 
taken to include minority populations of hot stars. 

5.1. Comments on other SPS models 

In this section we comment on a number of other SPS mod- 
els in the literature. The main purpose of the following discus- 
sion is to highlight the fact that many SPS models rely on the 
same ingredients as the models considered in detail herein, 



and therefore conclusions reached in earlier sections should 
be applicable to other models as well. 

The most p opular SPS mo d els no t considered herein are 
Starburst99 dLeifherer et al.l Il999t IVazauez & Leifhereri 
200 5]) and P eg ase dFioc & Rocca-Volmerangd 119971 : 
Le Bo rgne et al.l 120041) . In its most current version, 
Starburst99 combines the Geneva models for massive 
stars and the Padova isochrones for intermediate and low 
mass stars. Starburst99 employs the TP-AGB models of 
Vassil iadis & Woodl (119931) . which are the same models used 
in BC03. The BaSeL stellar library is used to convert the 
models int o observable SED s. For hot stars, the theoretical 
models of Smit h et ail d2002l) are used. 

For stellar populations not dominated by massive stars, 
Starburst99 will produce results similar to BC03. The treat- 
ment of nebular emission and massive stars, including Wolf- 
Rayet stars, is more sophisticated than either FSPS, M05, 
or BC03, and should therefore continue to be the model of 
choice when emission from massive star s is important. The 
recent Popstar model ( Molla et alj|2009l) also contains a so- 
phisticated treatment of massive stars, including the theo- 
retical atmosphere models for O and Wolf-Rayet stars of 
Smith et al.1 (120021) . and will therefore produce results simi- 
lar to Starburst99. 

Pegase utilizes the Padova isochrones and couples these to 
two different stellar spectra l libraries: BaSeL, and ELODIE 
(Prugniel & Soubiran 2 00l. ELODI E is an empirical library 
of 1500 stars. In lLe Borgne et al.l (120041) . various spectral 
indices are computed with Pegase and compared to observa- 
tions. The comparisons focus on a small sample of indices (in- 
cluding H/3, Fe5270, Mgb, and Fe5335) for Elliptical galaxies 
and MW star clusters. It would be fruitful to consider more 
comprehensive calibrating comparisons for the Pegase model, 
such as is considered herein. In any event, the underlying 
isochrones are similar to those used in BC03, so the predic- 
tions of Pegase should be similar to BC03. 

The Vazdek is mod els are described in Vazd ekisl dl999l) . 
IVazdekis et al.l d2003l) . and Vazdekis et al. (in prep). They 
utilize the Padova isochrones coupled to the empirical Miles 
stellar library. For old populations this model should thus be 
essentially identical to the FSPS predictions when the unmod- 
ified Padova isochrones are used. The Vazdekis model may 
thus fail to reproduce the observed near-IR colors of MC star 
clusters, since, as we demonstrated, the unmodified Padova 
isochrones are a poor match to those data. 

Significant efforts are currently underway to generate 
SPS models that handle a-enhancement self-consistently, 
in the sense that both the stellar evolution calculations and 
the a tmosphere models h ave the same ab undance patterns 
(e.g., ICoelho et al.l 120071: iLee et all l2009d) . These models 
will provide valuable constraints on the relative effects of 
a-enhancement. However, most models currently only fol- 
low stellar evolution through the end of the RGB, and so 
they cannot be used on their own to build mo del galaxies 
(the exception to this is the BaSTI database; Cordi er et alj 
I2007t iPercival et ail 120091) . We anticipate that the treatment 
of a-enhancement as a variable on par with age and Z in SPS 
models will likely become commonplace in the near future. 

From the above discussion it should be clear that the major- 
ity of SPS codes use the Padova isochrones with the BaSeL 
stellar library. Our comparisons and evaluations in previous 
sections, which rely on these ingredients, should therefore be 
applicable to many other SPS models. 
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6. SUMMARY 

We now summarize our principle results. 

• The latest stellar models from the Padova group can- 
not fit the near-IR photometry and surface brightness 
fluctuations of star clusters in the MCs. With our flex- 
ible SPS (FSPS) code we have modified the TP-AGB 
phase in the Padova isochrones to produce better fits 
to the data. The M05 models also fail to reproduce 
the data, both because the colors are too red and the 
age-dependence is incorrect. The BC03 and the FSPS 
model using the BaSTI isochrones fare well, although 
both models are somewhat too blue in the near-IR. A 
significant limiting factor in more accurate calibrations 
is the difficulty in constructing reliable photometry of 
MC clusters, owing to the fact that TP-AGB stars are 
luminous and rare. 

• At low metallicities (log(Z/Z Q ) < -0.5) and old ages, 
all models predict similar UBVRIJHK colors (the typi- 
cal variation between models is ks 0.05 mags). When 
compared to MW and M31 star clusters, the models 
are 0.2 mags too blue in V — K at log(Z/Z Q ) > -1.0, 
and 0.1 mags too red in J- K at all Z. The origin of 
the discrepancies in the near-IR is unclear. Decreas- 
ing the temperature of the RGB by ~ 200K alleviates 
much of the disagreement. If this modification is iso- 
lated to the brightest stars, then such a modification 
could improve agreement with integrated photometry 
without compromising agreement in CMD-space. The 
FSPS model with the BaSTI isochrones performs worst 
at low metallicities and old ages of the models consid- 
ered (although the disagreement is not dramatic), while 
the FSPS model with the Padova isochrones performs 
well. 

Multi-color CMDs of two metal-rich clusters, NGC 
6791 and NGC 188, are well-fit by FSPS. Spectral in- 
dices of MW star clusters are generally well-fit by both 
FSPS and BC03, although the models predict D„4000 
strengths too large and H5a strengths too weak com- 
pared to two log(Z/Z ) ta 0.0 clusters. 

• The ultraviolet photometry of the MW, M87, and M31 
star clusters can be well-fit by FSPS because FSPS con- 
tains flexible treatments of the post-AGB and horizon- 
tal branch evolutionary phases. This flexibility is essen- 
tial given the substantial scatter in the ultraviolet data 
at fixed metallicity, and given our inadequate theoreti- 
cal understanding of these phases. The M05 and BC03 
models perform less well because of their incomplete 
treatment of these advanced evolutionary phases. 

• SPS models are compared to ugrzYJHK photometry 
of massive red sequence galaxies. The FSPS, BC03, 
and M05 models fare well in most respects, with a few 
exceptions: all models are too blue in J — H; the M05 
model is somewhat too blue in Y—K; all models are far 
too red in g - r and u - g. These conclusions hold un- 
der the assumption that red sequence galaxies are com- 
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posed exclusively of old metal-rich stars with canonical 
evolutionary phases. The disagreement in the ugr col- 
ors can be alleviated if some combination of young (~ 1 
Gyr) stars, metal-poor stars, and blue straggler stars are 
added in small amounts. 

• Optical spectral indices of massive galaxies are gener- 
ally well-fit by the FSPS and BC03 models, with the 
important exceptions that both models predict D„4000 
in excess of observations and neither fit the observed 
locus of galaxies in the 004000-11(5/1 plane. These con- 
clusions hold for the same assumptions noted above, 
and, as above, small additions of some combination of 
young stars, metal-poor stars, blue straggler and hori- 
zontal branch stars can remedy this disagreement. 

• The FSPS and BC03 models adequately describe the 
grizYK colors of K+A, or 'post-starburst' galaxies, 
while the M05 model performs less well. Such galaxies 
contain a large proportion of intermediate age (0.5 < 
t < 2 Gyr) stars and thus provide unique constraints on 
the importance of TP-AGB stars in galaxies. 
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